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Recently, a self-contained trajectory-based formulation of non-relativistic quantum mechanics was 
developed [Ann. Phys. 315, 505 (2005); Chem. Phys. 370, 4 (2010); J. Chem. Phys. 136, 031102 
(2012)], that makes no use of wavefunctions or complex amplitudes of any kind. Quantum states are 
represented as ensembles of real- valued quantum trajectories that extremize a suitable action. Here, 
the trajectory-based approach is developed into a viable, generally covariant, relativistic quantum 
theory for single (spin-zero, massive) particles. Central to this development is the introduction 
of a new notion of global simultaneity for accelerated particles — together with basic postulates 
concerning probability conservation and causality. The latter postulate is found to be violated by 
the Klein-Gordon equation, leading to its well-known problems as a single-particle theory. Various 
examples are considered, including the time evolution of a relativistic Gaussian wavepacket. 



I. INTRODUCTION 

In this document, we develop a new formulation of 
single-particle relativistic quantum mechanics. Tradi- 
tionally, the formulation of quantum mechanics proceeds 
via a set of postulates, 1-5 which we do not find it neces- 
sary to repeat here. We do note, however, that the order, 
precise content, and even total number of quantum pos- 
tulates, vary from one treatment to the next. This situ- 
ation might be taken as an indication of the controversy 
or uncertainty that still exists — particularly around those 
postulates having to do with quantum measurement. The 
latter are especially nettlesome in the context of rela- 
tivistic quantum mechanics — where, e.g., it may not be 
entirely clear how to reconcile the "instantaneous" col- 
lapse of the wavefunction with subluminal propagation. 
On the other hand, the traditional quantum treatments 
all do agree on the first and most important postulate — 
that the state of a system be completely described by the 
quantum wavefunction, 'P. 

The wavefunction has thus always enjoyed a hallowed 
status in quantum mechanical theories, despite much his- 
torical and ongoing disagreement about its precise in- 
terpretation or physical significance. 6-19 Even those "al- 
ternative" interpretations of quantum mechanics that 
dare to challenge the first postulate still tend to respect 
the supremacy of "P. An example would be Bohmian 
mechanics, 7-13 which adopts a hybrid ontology wherein it 
is the wavefunction plus a quantum trajectory, together, 
that are needed to completely specify the quantum state. 
For many physicists, it may be difficult to even conceive 
of a quantum theory that makes no direct or indirect re- 
course to wavefunctions — after all, appears in every 
one of the five (or six) standard postulates. Neverthe- 
less, exactly just such a theory was recently formulated 
for non-relativistic quantum mechanics. 20-26 

For a number of reasons, it makes sense to try to extend 
the previous work to the relativistic case. As presented 
in this document, this goal is now also achieved — at least 
in the context of a single, spin-zero, massive, relativis- 
tic quantum particle, propagating on a flat Minkowski 



spacetime, with no external fields. Such a system might 
represent, e.g., a single Higgs boson particle — apropos 
to which, the recent news from CERN 27 is serving to 
stimulate demand for new and fresh approaches. In 
any case, generalizations for curved spacetimes, exter- 
nal fields, photons, particles with spin, multiple parti- 
cles, etc., together with detailed analyses of symmetry 
and stability properties, are planned for the future (al- 
though for each such development, a varying degree of 
required effort is anticipated). 

To be clear, the present approach is not a quan- 
tum field theory (QFT), but rather, a conserved par- 
ticle approach — in that sense, like the non-relativistic 
time-dependent Schrodinger equation (TDSE). As a fun- 
damental theory of relativistic quantum mechanics, it 
is safe to say that a particle-based strategy has been 
largely abandoned for many decades (with some no- 
table exceptions 28-30 ). Of course, the reasons for this 
date back to the earliest attempts to "relativize" the 
TDSE, starting with the Klein-Gordon (KG) equation 
in 1926. 3 < 5 < 10 > 31 - 37 Whereas the TDSE is first order in 
time and second order in space, the KG equation (which 
treats space and time on an equal footing) is second or- 
der in both — wherein lies the source of most of its diffi- 
culties. In particular, this leads to: (1) negative energy 
solutions, as well as (2) negative (indefinite) probability 
densities. 3 - 5 ' 10 - 34-37 

In 1928, Dirac improved matters somewhat, with the 
development of his famous first order (in both time and 
space) but multi-component equation, describing spin 
1/2 particles. 5 ' 10,38 ' 39 Dirac effectively solved problem 
(2), but not problem (1). By 1934 however, the "real" so- 
lution to both problems was hit upon — i.e., second quan- 
tization, and the development of QFT. 5,40-42 Though ob- 
viously serving us well in the ensuing decades, particu- 
larly for processes involving the creation/annihilation of 
particles, it can be argued that QFT introduces as many 
problems as it solves (e.g., pertaining to causality 41,42 as 
well as renormalization) , and in any event greatly compli- 
cates matters, both theoretically and conceptually. Pre- 
sumably, a viable, rigorous, single-particle theory of rel- 
ativistic quantum mechanics would therefore still be wel- 
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corned with open arms. 

The time is now ripe to revisit this notion. What 
makes us imagine that we can succeed where Klein, Gor- 
don, and Dirac evidently failed? The crucial development 
is the recent wavefunction-free reformulation of non- 
relativistic quantum mechanics, alluded to above. 20-26 
This approach is trajectory based, and in that sense remi- 
niscent of Bohmian mechanics. Unlike the Bohm theory, 
however, here, the traditional wavefunction, ^f(t, x), is 
entirely done away with, in favor of the trajectory en- 
semble, x(i, C) (where C labels individual trajectories) as 
the fundamental representation of a quantum state. The 
ensemble is such that there is exactly one trajectory pass- 
ing through every point in space (x) at any given point 
in time (t). The ensemble satisfies its own partial differ- 
ential equation (PDE) that is mathematically equivalent 
to the TDSE — even though it is nonlinear, real valued, 
second order in t and fourth order in C. It also satisfies 
an action principle, for which the Lagrangian consists 
of the usual classical part, plus a quantum contribution 
that incorporates intertrajectory interactions (i.e., partial 
derivatives of x with respect to C). 

From the perspective of developing a relativistic gen- 
eralization of quantum theory, the trajectory-based ap- 
proach is extremely well suited to the task. As in stan- 
dard relativity theory (and non-relativistic classical me- 
chanics), quantum solution trajectories are obtained as 
those that extremize the appropriate quantum action 
quantity- - whose relativistic form may be guessed as the 
Lorentz-invariant (or generally covariant) version of the 
non-relativistic quantum action described above. As in 
relativity theory, also, the quantum solution trajectory 
PDEs are inherently nonlinear. Although in the non- 
relativistic case, these happen to be equivalent to a linear 
wave PDE (i.e., the TDSE for there is no a priori rea- 
son to expect such a relationship to hold in the relativistic 
case. Indeed, if a viable (i.e., non-KG) relativistic linear 
wave PDE describing individual quantum particles were 
possible, then it probably would have been discovered 
decades ago by one of the aforementioned luminaries... 

At any rate, in this work we develop a relativistic 
generalization of our earlier non-relativistic trajectory- 
ensemble-based theory for quantum particles. It must be 
emphasized that new physics is being predicted here — 
which, in principle, could be validated or refuted by 
comparison with experiment. Because we are operat- 
ing largely in uncharted waters, it is possible that the 
present form of our equations may have to be modified 
(as was famously the case with Einstein's own equations); 
at the same time, however, general covariance considera- 
tions greatly restrict the form that such alternate dynam- 
ical laws might take. In any event, the form presented 
here is likely the simplest and most reasonable. We note 
that although only special relativity (SR) per se is con- 
sidered here, the mathematical development of our ap- 
proach relies heavily on curvilinear coordinate systems — 
and therefore, on the framework and tools of general rel- 
ativity (GR). 33,43 



Central to our approach is our (evidently) new defini- 
tion of simultaneity for accelerated particles. Even in SR, 
there is no good notion of the set of all spacetime events 
that occur simultaneously (from the particle's perspec- 
tive) with a given event on the particle's worldline, if it 
is accelerating. In this context, simultaneity is sometimes 
defined in the usual unaccelerated manner — i.e., using lo- 
cal inertial or comoving frames. 33 ' 43,44 This strategy fun- 
damentally fails, however, because it predicts multiple 
reoccurrences of the same spacetime events (e.g. pivot 
points 44,45 ), as well as the incorrect time ordering for 
distant, timelike separated events. Quite remarkably, the 
present, trajectory ensemble approach allows for a nat- 
ural and rigorous generalization of the simultaneity con- 
cept for an arbitrarily-moving particle — essentially, be- 
cause the quantum nature of the particle imparts a global 
character to it. 

The new relativistic quantum trajectory PDEs can be 
analyzed in various ways. In inertial coordinates, one ob- 
tains a form that is similar to the KG equation — yet dif- 
fers in one very crucial respect (discussed in Sec. VE). In 
retrospect, from a trajectory ensemble vantage, one can 
see clearly exactly where KG "got the physics wrong," 
in their attempt to relativize the TDSE. However, in or- 
der to do so, one must transform from the inertial frame 
to a certain curvilinear (albeit naturally arising) coordi- 
nate system, in terms of which the new PDEs are fourth 
order in space, and only second order in time (i.e., just 
like the non-relativistic quantum trajectory PDEs). The 
ramifications of this — and more generally, of the apparent 
inherent nonlinearity of the new PDEs — are not yet en- 
tirely clear. Thus, it may turn out that the present form 
is not always stable (or that other unanticipated prob- 
lems may eventually manifest) — although instability has 
not yet been observed, e.g., in the examples considered in 
Sec. VI. However, if such difficulties were to arise in the 
future, the author's view is that it should serve as a call 
to arms to look for a suitable rectifying modification of 
the present formulation — rather than as a condemnation 
of the general approach presented here, which seems to 
have much to offer. 

This document is quite long and comprehensive, as 
the requisite formal development is rather involved. We 
thus provide here a detailed overview of the remaining 
sections, with a specification of those subsections that 
might be skipped on an initial reading. Sec. II ad- 
dresses the basic mathematical structures that underpin 
the trajectory-based approach, in the relativistic context 
of a Lorentz-invariant four-dimensional (4d) spacetime. 
Sec. II A mainly establishes the notation and conventions 
as adopted here, but also introduces the trajectory en- 
semble four-velocity field; much of this material is stan- 
dard, and can be skipped by one versed in SR theory. The 
all-important simultaneity submanifold is then promptly 
constructed in Sec. II B. In Sec. II C , the "ensemble 
time" parameter is introduced, as a label for the different 
simultaneity submanifolds; this is found to be incompati- 
ble with the usual proper time, for reasons related to the 
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famous twin "paradox." 

Sec. Ill introduces the natural curvilinear coordinate 
system alluded to above (Sec. Ill A), together with var- 
ious probability density quantities. The most relevant 
equations in Sec. Ill A are Eqs. (12), (13), and (15), 
where the last two describe the form of the metric ten- 
sor in natural coordinates. Sec. IIIB introduces the first 
postulate of the trajectory-based approach (pertaining 
to probability conservation) and discusses the spatial 
scalar probability density, whereas Sec. Ill C considers 
the scalar invariant, 4d, and flux four- vector generaliza- 
tions. Sec. HID derives the covariant continuity equa- 
tion, and Sec. HIE discusses a particularly useful set of 
natural coordinates. The last three subsections of Sec. Ill 
are not as critical for an initial reading. 

Having laid out much of the mathematical framework 
in Sees. II and III, Sec. IV addresses dynamical consid- 
erations. The early part of Sec. IV A is critical, as it 
introduces the second postulate of the trajectory-based 
approach, pertaining to causality. This sensible condi- 
tion is satisfied by standard classical and non-relativistic 
quantum mechanics (TDSE), but not — it turns out — by 
the relativistic KG equation. Sec. IV B is a somewhat 
technical exposition on time reparametrization that may 
be skipped on an initial reading, whereas Sec. IV C is a 
review of the previous non-relativistic trajectory-based 
formulation, couched in the covariant language of GR. 

The meat of the new theory is presented in Sec. V. 
The new relativistic quantum PDE itself is readily ob- 
tained in Sec. VA [Eqs. (39) and (61)], although not in 
a form that is practically useful. That this PDE satisfies 
the principle of least action is established in Sec. VB, 
for those who wish to see how this comes about. The 
final, more practical form of the PDE [Eq. (77)] is then 
derived in Sec. VC — in which, also, an unexpected con- 
nection is established between the quantum and gravi- 
tational potentials. Sec. VD considers various limiting 
forms of the PDE (e.g., the non-relativistic limit), and 
Sec. VE presents a detailed comparison with KG the- 
ory; both may be skipped on an initial reading. Various 
examples are presented in Sec. VI, with the relativistic 
Gaussian wavepacket of Sec. VIC the most relevant. Fi- 
nally, a concluding discussion is provided in Sec. VII. 



and the metric tensor rj a p is the usual Minkowski one, 



V = 




(1) 



Note that we adopt the -H — | — h convention for the met- 
ric signature; also, factors of c are always explicitly in- 
dicated. Thus, the line element ds 2 = r] a p dx a dx /3 has 
units of length squared, whereas the proper time, r, de- 
fined via 



dr 2 = — f] a B dx a dx 13 , 
c z 



(2) 



has units of time. The Greek indices a, fi, 7, S, run 
over the four spacetime inertial coordinate labels, i.e. 
0,1,2,3, whereas fj,, v 1 k, a serve a similar function for gen- 
eral curvilinear coordinate systems (denoted X^). Latin 
indices run over spatial (or spacelikc) coordinate labels 
1,2,3, as per usual, with I, m (not to be confused with 
mass) used for inertial coordinates, and i, j, and k for 
general coordinates. 

A path ensemble (candidate solution trajectory ensem- 
ble) is uniquely specified via a contravariant vector field 
U a (in inertial coordinates), satisfying 



V, 



,p U a UP = -c 2 



(3) 



everywhere. Equation (3) above implies that U a is every- 
where timclike, and may be interpreted as a four-velocity 
field, i.e. 



U a = 



dx a 
Ik 



(4) 



Integration of Eq. (4) then gives rise to a set of Id sub- 
manifolds that foliate M, and therefore do not cross 
(even self-intersections are prohibited by the topology; 
the submanifolds are inextendible). This family of time- 
like curves thus constitutes the ensemble of paths, or can- 
didate solution trajectories. 

There is exactly one path passing through every point 
p in M; also, a one-to-one correspondence exists between 
paths (whose codimension is 3) and spatial points (the set 
x for fixed t). Moreover, all of the above properties are 
preserved under Lorentz transformations of the inertial 
frame, i.e. 



II. BASIC MATHEMATICAL STRUCTURE 



= A° 



(5) 



A. Preliminaries 

Let M be a 4d Reimannian manifold, representing the 
spacetime of a single, spin-zero, relativistic quantum par- 
ticle of mass m. For purposes of this study, M is pre- 
sumed flat (Ricci scalar R = 0). A global inertial frame 
can therefore be defined, in terms of which the inertial 
coordinates are x a — (ct,x) = (ct,x l ) — (ct,x,y, z), 



Specifically the transformed four-velocity field, U a , sat- 
isfies primed versions of Eqs. (3) and (4), and a cor- 
respondence can be established between paths and x' 
points, for fixed t'. 

We next introduce a set of path labeling parameters, 
C = C l , which uniquely identify individual paths within 
the ensemble. The C l are not yet coordinates per se, al- 
though later we will construct curvilinear coordinate sys- 
tems from them. Note that the C % values do not change 
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under Lorentz transformations. Since the labeling pa- 
rameters can (but in general will not) be identified with 
x at some specific t (for some specific inertial frame), we 
regard the C l as spacelike parameters. In any event, for a 
given path ensemble and choice of inertial frame, x(i, C) 
and C(i,x) are well-defined inverse diffeomorphisms of 
each other, for fixed values of t. 

Note that the above claims are subject to certain 
caveats, such as the possible existence of a measure-zero 
set of exception points. 22 

B. Simultaneity submanifolds 

Consider the tangent space T p for a point p in M. A 
3d spacelike "orthogonal subspace" can be defined as the 
set of all vectors W@ in T p that are orthogonal to U a : 

Vaf, U a W f3 = (6) 

The orthogonal subspace is a linear vector space in its 
own right; the collection of all such subspaces for ev- 
ery point p in M forms a subbundlc. We presume that 
the orthogonal subspaces can be integrated outward from 
the point p, to construct a corresponding integral sub- 
manifold. Specifically, this is a 3d embedded submani- 
fold of M, a spacelike hypersurface, that intersects ev- 
ery path exactly once. By starting this procedure from 
each point p that lies along a particular reference path, a 
one-parameter family of such hypcrsurfaces may be con- 
structed, which are non-intersecting, and otherwise foli- 
ate M. By construction, at every point p in M, the four- 
velocity vector U a is normal to the hypersurface passing 
through that point. 

We hereby refer to the above 3d hypersurfaces as "si- 
multaneity submanifolds." This terminology is justified 
through the following physical arguments. Consider a 
particle on a worldline passing through event p. At that 
instant, the velocity four-vector U a defines the local for- 
ward time direction for that particle. Likewise, the or- 
thogonal subspace defines the local spacelike directions 
for that particle — that is to say, the local set of events 
that occur simultaneously with p, from the particle's per- 
spective. This much, at least, is consistent with the idea 
of local inertial frames, and more importantly, with the 
Einstein Equivalence Principle. 33 ' 43 The problem in stan- 
dard relativity theory, of course, is that of extending 
these local notions of simultaneity in a global manner. 
Using inertial coordinate frames, this can be achieved 
for the special case of unaccelerated motion, but it is 
problematic for accelerating particles (a limitation that 
in hindsight, should perhaps seem a bit odd). 

In any event, the quantum relativistic theory de- 
veloped here provides a global concept of simultaneous 
events for a single, arbitrarily-moving particle, in the 
form of the simultaneity submanifolds described above. 
This is perhaps most physically meaningful if one adopts 
a "many worlds" -type ontological interpretation of the 
multiple particle paths/trajectories, according to which 



each trajectory worldline literally represents a different 
world, as has been discussed in previous work. 22 ' 24 ' 25 
The one particle is thus comprised of many "copies," dis- 
tributed across all space. Locally, the structure of the 
orthogonal subspace described above ensures that each 
particle copy agrees with its nearest neighbors as to which 
events occur simultaneously. Because of the global dis- 
tribution of copies, however, this notion can be extended 
globally throughout all of space. In this fashion — and 
much like Einstein's own orthogonal-ruler-and-clock con- 
struction of inertial frames — one builds a global manifold 
of simultaneous events, agreed upon by all particle copies, 
regardless of whether some or all of those copies are accel- 
erating, and despite the fact that they never cross paths. 

Note that the simultaneity submanifolds are not abso- 
lute, in the sense of being agreed upon by all observers. 
A different quantum observer (or particle) would have 
its own copies, its own trajectory ensemble, and (in gen- 
eral) an entirely different set of simultaneity submani- 
folds. This is as it should be. Note that for two quantum 
observers to agree completely on simultaneity, every tra- 
jectory in one ensemble must match the corresponding 
trajectory in the other ensemble (i.e., the two trajectory 
ensembles must be identical). 

Likewise, a "quantum inertial observer" is char- 
acterized via an ensemble of parallel straight-line 
trajectories — i.e., by 

U a = const". (7) 

Thus, for example, if x a satisfies Eq. (5), then the con- 
tours of x° (expressed in the x a coordinate system) de- 
fine the simultaneity submanifolds for a quantum iner- 
tial particle, whose trajectory orbits are given by the in- 
tersections of the x l contours. As indicated in Fig. 1, 
this special case is of course, entirely consistent with 
the usual global notion of simultaneity for unaccelerated 
particles — except that here, it is obtained in a more rig- 
orous, essentially completely local fashion. The reason is 
that U a is a four- velocity vector field — defined on all of 
M, rather than just along a single trajectory. 

More generally — i.e., for accelerated motion — the rel- 
ativistic quantum trajectories and simultaneity subman- 
ifolds behave more along the lines indicated in Fig. 2. 
Unlike comoving frames, 33,43 ' 44 the simultaneity subman- 
ifolds of the present theory are curved. This is appropri- 
ate, given that the trajectories are also curved (albeit 
in the extrinsic rather than intrinsic sense). The curva- 
ture of the simultaneity submanifolds also enables them 
to avoid intersecting each other — thus also avoiding the 
problems of multiple reoccurences and incorrect time or- 
derings that plague the comoving frame approach. Fi- 
nally, the simultaneity submanifolds are everywhere or- 
thogonal to the trajectories, and thus each submanifold 
consists only of spacelike separated events. 
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Figure 1 
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FIG. 1: Contour plot of Id inertial coordinates, x a = 
(ct' , x'), as functions of the inertial coordinates, x a = (ct, x), 
to which they are related via a Lorentz transformation 
[Eq. (5)]. The dashed lines are the contours of t 1 (ct, x). These 
represent "simultaneity submanifolds" — i.e., sets of spacetime 
events that occur simultaneously — for an inertial observer 
moving with the primed frame. The solid lines are the con- 
tours of x'(ct,x), any one of which could represent the tra- 
jectory of the moving inertial observer or particle, according 
to standard SR theory. Collectively, the solid lines constitute 
a relativistic quantum trajectory ensemble, describing a single 
relativistic particle undergoing quantum inertial motion. 



C. Ensemble time, and the generalized twin 
"paradox" 

It should be noted that the ability to construct global 
simultaneity submanifolds is not automatic, but in fact, 
induces a slight constraint on the allowed form of U a . 
Using Frobcnius' theorem, 33 ' 46 and the fact that the U a 
field is (presumed) smooth, it can be shown that only 
fields for which h U a is four-curl-free (for some scalar field 
h) are permissible. (Technically, this is slightly more re- 
strictive than the actual condition, but it is a sufficient 
condition that suits our purposes better). This condi- 
tion implies that the four-velocity field is the normalized 
gradient of some scalar field A: 



U a =c 



-if 



(8) 



d/dx a (and it is presumed 



where n" 13 — r] a p and d a 
that d X > 0). 

The quantity A is a timclike parameter that we call an 
"ensemble time." Actually, it is a full-fledged scalar field, 
and can therefore be interpreted as a global time coordi- 



nate. Eqs. (6) and (8) imply that the A = const contours 
are the simultaneity submanifolds. The term "ensemble 
time" is therefore justified, as all members of the ensem- 
ble (i.e., all particle copies) agree that events correspond- 
ing to the same value of A occur simultaneously. The 
actual A value itself, however, is not uniquely defined. In 
general, any transformation of the form A — > A' = A' (A) 
yields a new ensemble time coordinate with the same 
contours, and which otherwise also satisfies Eq. (8). For 
the moment, we treat all such choices equally. Later, we 
will identify special candidates to serve as the "ensem- 
ble proper time," A = T (Sees. HIE and VC) — i.e., the 
ensemble analog of the usual single-path proper time, r. 

Even for a path ensemble, it is straightforward — and 
often convenient — to construct a true proper time coor- 
dinate, r, as a scalar field on M. For each path in the 
ensemble, labeled by C, one simply chooses a reference 
point p at which r is taken to be zero, and then integrates 
Eq. (2) along the path to find the value of r at all other 
points. There is thus a freedom in the definition of the r 
coordinate, associated with the particular choice of ref- 
erence point p for each path. This freedom corresponds 
to coordinate transformations of the form 



r ->t' = t'(t,C) = r + A(C), 



(9) 



where the shift, A(C), varies across paths. 

Although the r coordinate is useful in its own right, 
note that any choice of r is in general incompatible with 
any A coordinate of the form described above. In par- 
ticular, for two different paths of a given ensemble, both 
starting at the same simultaneity submanifold (A con- 
tour) and ending at another, the elapsed ensemble time 
AA is the same, but the elapsed proper times At are 
(generally) different. Thus, r per se cannot generally 
serve as a good ensemble time coordinate, A. 

The situation above is not problematic, and in fact, 
reflects nothing more than a generalized version of the 
well-known (but poorly named) twin "paradox." The 
conventional twin paradox has one twin leaving the other 
at spacetime point p , only to rejoin him or her at a 
later point pf, after having undergone accelerated mo- 
tion. The accelerated trajectory of the second twin is 
necessary, in order that their paths may recross (the 
stay-at-home twin is presumed to undergo inertial mo- 
tion). Since the starting point for both twins is in fact 
the same event po, there can be no question but that the 
departures occur simultaneously. Likewise, the reunion 
at pf is the same event for both twins, and must there- 
fore also occur simultaneously. Nevertheless, we know 
that the elapsed proper times for the two twins are dif- 
ferent, with the stay-at-home twin having aged (in some 
renditions, very significantly) more than his or her more 
adventurous sibling. 

A similar situation characterizes our ensemble of 
paths — except that we now have a global definition of 
simultaneity, so that we can specify that two such paths 
begin "at the same moment" even if the initial spacetime 
events arc different. Likewise, we can uniquely specify 
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the simultaneous "end" of the two paths, as the points 
where these paths intersect a different A contour. The 
global property of the simultaneity submanifolds is in- 
deed required, as any two paths within a given path en- 
semble never cross. Even though the two paths start 
and end simultaneously, the generalized twin paradox im- 
plies that there is no reason to expect the two elapsed 
At values to be the same — and indeed they are not, 
in general. For the special case where all paths are 
moving inertially — i.e., for a "quantum inertial observer" 
(Fig. 1) — then the elapsed proper times for all paths are 
equivalent, and it is permissible in this instance to use 
r as an ensemble time coordinate, A. We thus make it 
a requirement of any reasonable definition for an ensem- 
ble proper time, T, that it should reduce to the usual 
proper time, r, in the special case of a quantum inertial 
observer. 



III. GENERAL COORDINATES AND 
PROBABILITY 

A. General (curvilinear) coordinates and natural 
coordinates 

All of the equations of Sec. II have been presented 
in a way that is manifestly covariant, with respect to 
arbitrary coordinate transformations (diffeomorphisms). 
The general (or curvilinear) coordinate version of all such 
results is obtained by replacing the indices a, f3, etc. with 
/i, v 1 etc., the inertial coordinates x a with general coor- 
dinates X^, and the Minkowski metric tensor with the 
generalized form, 

dx a Qx^ 

9^ = V<*0 q^j (10) 

Note also that in principle, all partial derivatives d a 
that appear in Sec. II should be replaced with covariant 
derivatives, 33 ' 43 denoted here as V M . However, these ap- 
pear only in Eq. (8) , where they are applied to a scalar in- 
variant field (A), for which it is well-known that V M = 
(in any coordinate system). 

It is sometimes convenient to write Eq. (10) in matrix 
form, 

9 = J T ■ V ■ J, (11) 

where J a ^ — d^x" is the Jacobian matrix for the co- 
ordinate transformation x a — > X^. We also define the 
inverse transform Jacobian matrix, K, as K^ a = d a X^ . 
Being true inverses of each other, J a ^K^p = 6p, and 
K^a J a v = 6£, where in this context, 5 is the Kronecker 
delta function. 

Of all of the general coordinate systems that could be 
used to characterize our flat Minkowski spacetime man- 
ifold, M, obviously the global x a inertial frame coordi- 
nates considered in Sec. II are a natural choice. How- 
ever, for a given path ensemble, other natural choices 



also arise, based on quantities that we have already in- 
troduced. Let us hereby define a system of natural coor- 
dinates as the curvilinear choice, 

X" = (cA,C) = {c\,&) = (c\,C\C 2 ,C 3 ), (12) 

where A and C % are defined in Sec. II. In that section, 
these quantities were considered parameters; however, it 
is clear from the discussion therein that they can be pro- 
moted to full-fledged coordinates, forming a good coor- 
dinate system under the conditions discussed. 

The physical meaning of the natural coordinates is 
such that the intersection of the contours of the three 
functions, C z (x a ), define the individual paths within the 
ensemble, whereas A itself describes the (ensemble) time 
evolution along a given path (see Figs. 2 and 3). Note 
that the introduction of a factor of c into the definition 
of X° is consistent with the interpretation of A as a time- 
like coordinate. Likewise, the C l coordinates, which (as 
we have seen) can be used to label points on a given 
simultaneity submanifold, can be regarded as spacelikc 
coordinates. That is not to say, however, that A must 
have units of time, and C % units of length. Rather, it 
is better to think of natural coordinates as being any 
choice that respects the natural division of M into space 
and time that is induced by a given path ensemble (as 
described in Sec. II). 

Further justification for this interpretation of the nat- 
ural coordinates is provided by the fact that for any 
coordinate transformation of a form that respects this 
time/space division — i.e., the reparametrizations A — > 
A' = A' (A) and C -> C = C'(C) — the new coordinates, 
X^ = (c A', C), are also seen to be natural coordinates, 
by the definition given above. For the moment, we treat 
all such choices equally. Later, however, after we have 
introduced suitable probabilistic structures and dynami- 
cal elements on M, we will find that a particular choice 
naturally emerges [i.e., X^ = (cT, P)]. 

An important feature of any set of natural 
coordinates — consistent with the above interpretation — 
is that the metric tensor g be block-diagonal. Thus, 
9oi = 9i0 = Vai3J a oJ f3 i = 0, and so 

#-(? ;)• o»> 

In Eq. (13) above, 7 — the so-called "spatial metric" — 
represents the 3x3 spatial block of the full metric tensor 
g. Natural coordinates are therefore orthogonal with re- 
spect to the division of space and time (although the 
spacelike coordinates C l are not necessarily orthogonal 
amongst themselves). This is a manifestation of the fact 
that the timelike and spacelike subspaces of the tangent 
vector spaces were constructed, by design, to be orthogo- 
nal to each other [Eq. (6)]. Note that reparametrizations 
of A do not affect 7, and reparametrizations of C do not 
affect goo- 

The block-diagonality assertion above is readily proven 
in terms of the inverse metric tensor, 
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FIG. 2: Quantum trajectories (solid curves) and simultaneity submanifolds (dashed curves) for the Id c~3 relativistic Gaussian 
wavepacket of Sec. VIC, as represented in the inertial coordinate frame, x a — (ct,x). At t = (i.e., along the x axis), 
all spacetime events occur simultaneously. The curvature of the simultaneity submanifolds at later times is an indication 
of relativistic dynamical effects, which are quite pronounced. The trajectories are everywhere normal (orthogonal) to the 
simultaneity submanifolds. The dark gray and light gray patches correspond to those in Fig. 3. 



g- 1 = K-f]-k T . (14) 

Note that K° a oc d a X, whereas the contravariant vectors 
r\ a ^K % p describe displacements within the simultaneity 
submanifolds. Thus, from Eqs. (6) and (14), we have 
g° l = g t0 = 0, and the same must therefore be true for g 
itself. Another useful set of relations to follow from the 
time/space orthogonality of the natural coordinates is 

--(£)'. < 15 > 

The last equality relates the time-time component of the 
(inverse) metric tensor to dX/dr, the local relation be- 
tween ensemble and proper time (along a given trajec- 
tory/path). 

Finally, we have a relation for the determinant of g, 
which will appear in various expressions. Following stan- 
dard convention, we take 

g = detg = -(dot J) 2 = -J 2 = -l/K 2 (16) 

No absolute value is implied in Eq. (16) above; indeed, 
we see that g must be negative. Equation (16) above 
is completely general. For natural coordinates, however, 
this reduces to 

9 = .9oo7. where 7 = detj (17) 



(and 7 should not be confused with the index). It is clear 
from the previous definitions that 500 < and 7 > 0; 
thus, the -+++ signature of the inertial frame metric is 
retained, in keeping with our time/space interpretation 
of the natural coordinates. 



B. Probability in 3d 

Classical statistical mechanics, trajectory-based quan- 
tum mechanics, and relativistic hydrodynamics all in- 
clude the notion of density fields that propagate through 
time and S p aC e. 10 < 1:L > 22 ' 25 > 26 ' 33 > 43 The densities are asso- 
ciated with corresponding flux vectors, that govern the 
local motion of fluid elements. Densities and their corre- 
sponding fluxes obey a continuity relationship, expressing 
the physical conservation law for a given quantity — be it 
energy/mass, charge, probability, or number of particles. 

As we seek here to generalize the non-relativistic the- 
ory of quantum trajectory ensembles, and since the non- 
relativistic theory is mathematically equivalent to the 
TDSE, the conserved quantity in question can be taken 
to be probability. In the non-relativistic case, estab- 
lishing the appropriate continuity equation is relatively 
straightforward — essentially because the time coordinate 
t is uniquely determined. In the relativistic case, how- 
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FIG. 3: Quantum trajectories (solid curves) and simultane- 
ity submanifolds (dashed curves) for the Id c~3 relativistic 
Gaussian wavepacket of Sec. VI C, as represented in a natural 
coordinate frame, X M = (c A, C). All velocity and flux vectors 
point vertically (e.g., J7 M ), and all force vectors point horizon- 
tally (e.g., fg). The probability-length contained in each gray 
patch is identical to that of the corresponding patch in Fig. 2, 
and independent of the choice of A coordinate. This value 
is not the same for the light gray and dark gray patches, 
although these two patches do have the same value for the 
spatial scalar probability density, /(C). 



ever, there are myriad probability-related quantities that 
can in principle be defined. One can construct both scalar 
and vector densities (current or flux four- vectors), as well 
as stress-energy-type tensors — all of which can be further 
subcategorized as to whether they exist on the full 4d 
spacetime, or only on submanifolds. In addition, whereas 
density quantities are generally not true invariants (i.e., 
with respect to general coordinate transformations), 33,43 
one can often construct invariant versions of these quanti- 
ties. To cut through the morass of possibilities, wc apply 
the standard procedure of contemplating how such quan- 
tities, and the relationships among them, should trans- 
form under various coordinate transformations. 

We also rely on the following assumption, borrowed 
from the non-relativistic theory: 

• Postulate 1: Probability is conserved along quan- 
tum trajectories. 

In the non-relativistic case, "probability" means the 
true, unitless, probability element — i.e., the (spatial) 
probability density times the (spatial) volume element, 
p(t, x(£)) d 3 x{t). Postulate 1 stipulates that along a given 
trajectory (i.e., for fixed C), this quantity is independent 
of t. Under the coordinate transformation x — > C (for 



fixed t), the probability density transforms as 
p(t,x(t))d 3 x = f(t,C)d 3 C, 



(18) 



because the probability element itself must be a scalar 
invariant. Since C and d 3 C are independent of t along a 
trajectory, the probability conservation postulate there- 
fore implies that f(t, C) = /(C) is itself independent of 
t. Thus, /(C) d 3 C assigns a definite probability element 
to each trajectory in the ensemble, for all time. 

We posit a similar situation in the relativistic case. 
Let /(C) be introduced as a probability density on C 
space — i.e., on the set of individual paths within a given 
ensemble. This density has units of C~ 3 , and is presumed 
to be normalized to unity: 



Jf(C)d 3 C 



1 



(19) 



The value of /(C) must be independent of the value 
of A, in order to satisfy the probability conservation 
postulate — but equally importantly, it must be invari- 
ant of A reparametrizations, A — > A' (A). The function 
/(C) can thus be "pulled back" to the individual simul- 
taneity submanifolds on M, 33,43 but in no sense should 
it be regarded as a probability density on M itself. Wc 
henceforth refer to /(C) as the "spatial scalar probability 
density" on C. 

The situation is exemplified by Fig. 3, which depicts 
the spacetime manifold M as charted in a natural coor- 
dinate system. In this frame, individual paths appear as 
vertical lines, and simultaneity submanifolds (surfaces of 
constant A) as horizontal lines. The evolution of ensemble 
time — technically, the A-parametrized family of M — > M 
diffcomorphisms generated by the A velocity vector field, 



dXV_ 
~dX~ 



(c, 0,0,0) 



(20) 



— serves to advance the simultaneity submanifolds ver- 
tically through the spacetime M. Of course, the "rate" 
at which this occurs depends on the choice of ensem- 
ble time coordinate A — which in turn affects the vertical 
length scale in Fig. 3, and the "thickness", cdX, of a four- 
volume element, cd\d 3 C. The changes in these quanti- 
ties brought about by a reparametrization of A would in 
turn affect the values of any 4d density quantities on M. 
However, A reparametrization has no effect on the simul- 
taneity submanifolds themselves, nor on the 3d quanti- 
ties that live on these submanifolds, such as the spatial 
scalar probability density, /(C), and the spatial volume 
element, d 3 C. 

On the other hand, we of course demand that 
/(C) should transform as a probability density under 
reparametrization of the spatial coordinates, C — > C: 



f'(C')d 3 C = f(C)d 3 C ; /'(C) 



-/(C) (21) 



Technically speaking, as a function of C, /(C) is not a 
scalar invariant field, but a 3d scalar density of weight 
W = -l. 
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C. Probability and flux in 4d 

As indicated, the 3d (spatial) /(C) density transforms 
in a well-prescribed way, under all coordinate transforms 
that preserve the natural coordinate structure. Our ul- 
timate goal, however, is a fully covariant formalism, for 
which any choice of coordinates might in principle be 
used. We should therefore also consider the appropriate 
4d scalar probability density (and the corresponding flux 
four- vector) . A natural way to achieve this is through 
the use of the invariant form of the spatial scalar prob- 
ability density — which is a true scalar invariant (weight 
W=0). We denote such invariant forms of density quan- 
tities with an asterisk superscript. The "spatial scalar 
invariant probability density" is thus found to be 



/*(C) = 



/(C) 

V7 ^ 



(22) 



as is evident from Eq. (21). 

The advantage of working with true scalar invariants 
is that it is quite straightforward to pull back a scalar 
invariant function on M to any submanifold of M, sim- 
ply by restricting the domain of the function accordingly. 
In the present case, we would like to define a 4d scalar 
invariant probability density on all M, denoted p*(X tl ). 
What we have is a 3d spatial scalar invariant probability 
density /* defined on individual simultaneity submani- 
folds; but since the set of all such submanifolds foliate 
M, /* can be promoted to an actual scalar invariant 
function on M, i.e. f*(X tl ). It is therefore completely 
natural to define 



p*(X») = f*(X») 



(23) 



in a generally covariant manner. The desired 4d scalar 
probability density (of weight W = — 1) then becomes 



P (xn = vw\p*(xn 



f[C(X»)] (24) 



Equation (24) above is the completely general covari- 
ant expression. In this context, however, it must be un- 
derstood that g refers to the completely general coor- 
dinates X^, whereas 7 refers to the coordinates C. In 
the special case where X^ = (c A, C) is taken to be the 
natural coordinate frame itself, then Eq. (24) reduces to 



p(X") = V^~o /(C) 



/(C). (25) 



In any event, p(X^) behaves as a true 4d scalar density 
should; the element p(X^) d A X is a scalar invariant. This 
element is not unitless, however; it has units of length, 
as the invariant form p* has units of length" 3 . 

From the scalar probability density p{X^), we wish to 
construct a flux four-vector, , such that 



1 1 
-c p . 



(26) 



Intuition (e.g., from electromagnetic current) suggests 
that should transform as a contravariant four-vector, 
which would make p a scalar invariant. In SR, this is com- 
pletely correct, in the sense that p'(x a ) = p(x a ) is unal- 
tered under Lorentz transformations, because the volume 
element d 4 x' = d 4 x is unchanged. Thus, indeed trans- 
forms as a four- vector under Eq. (5). In the completely 
general context however, incorporating transformations 
to curvilinear coordinate frames, one finds that p and 
should transform, respectively, as scalar and contravari- 
ant vector densities, of weight W=—l. This distinction 
is quite critical from the perspective of the covariant con- 
tinuity equation (to be defined shortly). 

From Eqs. (3) and (26), and the requisite transforma- 
tion properties, the desired flux quantity is evidently 



f = P (xn u» - /(c) w, 



(27) 



where the first equality is generally true, and the second 
holds [from Eqs. (20) and (25)] in a natural coordinate 
frame. 



D. Covariant continuity 

Next, we derive a continuity equation in a natural co- 
ordinate frame. We expect this to take the form of a 
vanishing four-divergence of j M , i.e. 



= 



(28) 



From Eqs. (20) and (27), the flux in a natural coordi- 
nate frame is given by j M = (c/(C), 0,0,0), for which 
Eq. (28) is clearly satisfied, by virtue of the probability 
conservation postulate. 

Is Eq. (28) also true for general coordinates? Since 
derivatives of vector quantities are involved, ordinar- 
ily, one would have to replace the partial derivatives in 
Eq. (28) with covariant derivatives, in order to obtain a 
generally covariant expression. However, for the special 
case of a contravariant vector density with weight W = — 1 
(or a covariant vector density with weight W = +1) 
it turns out that the covariant four-divergence is equal 
to the usual (partial derivative) four-divergence. Thus, 
c>„ j M = holds in all coordinate frames, and there- 

fore, so does Eq. (28). We take this to be the covariant 
continuity equation. 

The four-divergence property described above can be 
proven as follows. Let j 11 be a contravariant vector den- 
sity of weight W = — 1. The standard expression for the 
covariant derivative is 



(29) 



where the arc the famous Christoffel symbols, 33 ' 43 
defining the (metric compatible) connection: 



r^ K = ^ a {d v g Ka + d K g av - d a g VK ) (30) 
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To obtain the covariant divergence from Eq. (29), v is 
set equal to /i, which then becomes a dummy index that 
is summed over (Einstein notation). By the symmetry 
of the Christoffcl symbols with respect to their lower in- 
dices, the second and third terms in Eq. (29) then cancel, 
leaving the desired equality. 

The covariant continuity equation is easiest to inter- 
pret in a natural coordinate frame. Here, the j M vector 
field is "vertical" (Fig. 3), and simply describes how prob- 
ability density is carried along under the action of the V 11 
velocity field. It is important to note that it is the spatial 
scalar probability density, /(C), that is conserved in this 
sense, not the scalar probability density p(X^) per se. 
To within a constant factor of c, the former is just the 
zeroth component of the flux four- vector, j°, in a natural 
coordinate frame. More generally, i.e. in an arbitrary 
coordinate system, it is (j°/c) that we ordinarily think 
of as the "probability density". In particular, if is 
any coordinate system for which X° and X % are orthogo- 
nal, then Eq. (28) and the generalized Stokes theorem 33 
imply that 



/ 



fiX") dX 1 dX 2 dX 3 = const, for all X° (31) 



[provided i Q (X^) vanishes asymptotically]. 

Note that the orthogonal condition above is not the 
same as the natural coordinate condition, and is in fact 
far less restrictive. For instance (and in addition to any 
natural coordinate frame) it includes any inertial coordi- 
nate frame, X^ = x a , in which case Eqs. (28) and (31) 
reduce to the usual SR forms. The inertial frame forms 
of the density and flux quantities bear further discussion. 
From Eq. (24), we find that 



p(x a )=p*(x a ) = 



f[C(x a ) 
\/7 



(32) 



Likewise, from Eqs. (25) and (27), we obtain 

r 



f[c] na „ f[C] 

u — c — , j o 



y/l ~ ~ V-.9oo7 
„/[C] dt 



(33) 



Clearly, p{x a ) and (j a /c) both have units of length -3 , 
as they should. Note that goo, 7, and J a o refer to the 
natural coordinate frame. 



E. Uniformizing natural coordinates 

Let us return to a natural coordinate frame, and con- 
sider the ensemble time evolution, as indicated in Figs. 2 
and 3. Although the probability element p(X^) d 4 X is 
not conserved along trajectories, it is invariant. Thus, 
the probability-length contained within the dark gray 
patches in Fig. 2 and Fig. 3 is the same, and the 
probability-length within the light gray patches in Fig. 2 



and Fig. 3 is the same, but the value for the dark gray 
patches is not the same as for the light gray patches. 
Fundamentally, this is because A and r are generally in- 
compatible, and so dr/dX must vary across M. Put an- 
other way, a reparametrization of A could expand the 
vertical scale of the light gray patches without necessar- 
ily changing that of the dark gray patches, which has a 
corresponding effect on the p values. Clearly, p cannot 
be conserved for a general choice of ensemble time, A. 

Nevertheless, it would be nice if we could define a 
"proper-time like" reparametrization of A, for which 
p(X^) was as close to being conserved as possible. The 
best we can manage along these lines is to choose a A 
such that the spatial integral of p(A M ), i.e. 



J P (xnd 3 c = J 



dr 



p(X»)d 3 C = I — f(C)d 3 C 



(34) 



is conserved over time. We can therefore define, as an 
ensemble proper time, T, that choice of ensemble time A 
for which 



/ 



dr_ 
df 



/(C) d 3 C = 1, for all T. 



(35) 



The interpretation of Eq. (35) is straightforward. By 
multiplying both sides by dT, we see that, at any time 
T, the interval dT is just the path ensemble average of 
the proper time interval dr. If the dr intervals happen 
to be the same for all paths in the ensemble at a given 
T, then dT = dr at that time. If this is true at all times, 
then T = t — as describes the special case of quantum 
inertial observers. In general though, the dr values dif- 
fer across the ensemble, and so dT is obtained as the 
probability-weighted average of the dr values for the in- 
dividual paths. Note that for a given path ensemble, 
the T coordinate as specified here is uniquely defined, to 
within an overall additive constant. 

As reasonable as the above prescription may appear, 
Eq. (35) is not the most natural choice for an en- 
semble proper time, T; that choice will be introduced 
in Sec. VC. For one thing, even with Eq. (35), the 
probability-length of the dark gray patches in Figs. 2 
and 3 is still not equal to that of the light gray patches — 
although the spatial integrals across the simultaneity sub- 
manifolds are now equal. The ensemble proper time as 
defined above can be regarded as a "uniformizing coor- 
dinate," in the sense that it does the best possible job 
of spreading out the scalar probability density, p(X^), 
uniformly throughout (ensemble) time. 

One might also consider a similar reparametrization of 
the (natural) spatial coordinates, C. Here, it is possible 
to perfectly uniformize the probability distribution /(C), 
and in fact, exactly this procedure has been used in the 
previous non-relativistic work. 25 We hereby denote such a 
uniformizing choice of spatial coordinates for C as C = P, 
defined such that 



/(P) = 1- 



(36) 
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For a single spatial dimension, the P coordinate is 
uniquely determined, apart from an additive constant 
(which fixes the range of allowed P values); it can be 
interpreted as the total (integrated) probability that ex- 
ists to the left of a given path in the ensemble (hence 
the nomenclature "P," often used in such contexts). 
In higher dimensional spaces, the P coordinates are 
also uniquely determined — apart from (spatial) volume- 
preserving diffeomorphisms. 

With T taken to be any ensemble proper time [e.g., 
either that of Eq. (35) or of Sec. VC], we regard = 
(cT, P) to comprise a set of "uniformizing natural coor- 
dinates," in terms of which T has units of time, P is unit- 
less, /(P) = 1 is unitless, and p(X») = (dr/dT)//^ 
has units of length" 3 (as for the inertial frame case). 



IV. DYNAMICAL CONSIDERATIONS 

We now have many of the elements in place that we 
need to develop a trajectory-based theory of relativis- 
tic quantum dynamics. However, there are still a few 
more issues, dynamical in nature, that must first be ad- 
dressed. In the first two subsections below, we consider 
single-trajectory classical dynamics, both non-relativistic 
and relativistic. Even in this context, there are certain 
subtle aspects that turn out to have extremely impor- 
tant ramifications for the relativistic quantum case. In 
the third subsection below, we consider the trajectory- 
based non-relativistic quantum theory, presenting some 
key results from earlier work, and setting the stage for 
addressing quantum effects in the relativistic context. 



A. Classical dynamical considerations 

Consider a single, non-relativistic classical particle. 
The dynamics are described by the trajectory, x (t). The 
velocity vector is x l = dx l /dt, and the non-relativistic 
classical equation of motion is 



mxl = JVit^ = fl 



dx l 



(37) 



where V(t,x) is the potential energy field, and f is 
the force vector. Note that we are considering the gen- 
eral, nonconservative case where V may depend on t as 
well as x. The reason is that we wish to adopt a 4d 
spacetime viewpoint, even in this non-relativistic con- 
text. We therefore continue to use notation such as u x a " 
and u g a p" mostly without ambiguity. Thus, for exam- 
ple, the Euclidean metric, g a p — 5 a p is presumed. 

From this 4d vantage, a striking feature of Eq. (37) is 
that the force vector does not include a time component, 
even when V depends on t. In other words, the dynamical 
equations do not make use of what would be the entire 
force four- vector, f a = f a = d a V, but only the spatial 
components of this vector — i.e., the projection onto the 



(non-relativistic) simultaneity submanifold (just 3d space 
itself, x). If the time component were used, it would lead 
to an additional dynamical equation of the form 



mi = f° oc - 



dV 

~dt' 



(38) 



which is manifestly incorrect for time-dependent poten- 
tials, because i = di/dt = d(l)/dt = 0. 

From a physical standpoint, the restriction to just the 
spatial components of the force vector is not surprising. 
Suppose that the partial time derivative of V(t, x) were 
employed in the determination of the dynamical force. 
That would mean that this force, which drives the fu- 
ture time evolution of the particle, would itself depend 
on the future states of that particle — a highly untenable 
situation vis-a-vis causality, which can be reasonably dis- 
missed. Of course, time reversibility implies that the 
force must also be independent of the past state of the 
particle, which leaves only the present. In this manner, 
we are led to: 

• Postulate 2: All force vectors, together with the 
quantities used in their construction, must "live" 
on the simultaneity submanifolds. 

The meaning of this statement will be made more precise 
as we go along. In any case, we take Postulate 2 as 
a necessary condition for any viable physical theory — 
at least for particles with mass, in the context of non- 
relativistic or relativistic classical or quantum mechanics. 
As shown above, it certainly holds for non-relativistic 
classical particles. 

It also holds for relativistic classical particles, as we 
now demonstrate. Of course, it is not possible to define 
a global simultaneity submanifold from a single relativis- 
tic trajectory, x a (r). All that is required for the present 
purpose, though, is a local specification of simultaneity, 
which we do in fact have. As discussed in Sec. II, this is 
found — for a given point p — in the orthogonal subspace 
of T p that is orthogonal to U a at p. According to Postu- 
late 2, the relativistic force vector /" must belong to this 
orthogonal subspace. That this is the case can be shown 
using the relativistic classical equation of motion: 



m 



dr 2 



ra- 



dii" 



r 



(39) 



(note that we are still working in an inertial frame). From 
Eqs. (3) and (39), 



to d(n a pU a U f3 ) 
~2 dr 



= V*p U a f = 0. (40) 



Thus, indeed, f a belongs to the orthogonal subspace, and 
hence points along the simultaneity submanifold. 

Let us pause to consider the ramifications of Eq. (40). 
Since the orthogonal subspace at a given point p depends 
on the trajectory, this implies that f& itself must also de- 
pend on the trajectory at p — i.e., not just on the point 
itself. Thus, for example, a relativistic force vector can 
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never be simply the gradient of some scalar invariant po- 
tential: f a 7^ —rj al3 dp$>. On the other hand, one simple 
way that Eq. (40) can be achieved is by taking f a to be 
the product of an antisymmetric tensor and U a . This 
is the general form for the relativistic force that ensues 
when the Lagrangian includes a term that is the scalar 
product of U a and some four-vector field. The classic ex- 
ample is the vector potential, A a , of electromagnetic the- 
ory, for which the resultant antisymmetric tensor, F a p, 
is the electromagnetic field strength tensor. Thus, the 
"velocity-dependent" nature of a relativistic force field — 
well-known in the electromagnetic context — is seen to be 
a general requirement arising from Postulate 2. 

Note that, whereas f a does indeed belong to the or- 
thogonal subspace of T p , the same cannot be said for A a 
and U a — despite the fact that these are clearly "elements 
used in the construction" of f a . Postulate 2 can still be 
regarded to be satisfied — albeit in a weaker sense than 
for f a itself — because the A a and U a vectors belong to 
the tangent space T p for the point p itself, and not some 
other point. Thus, in that sense, these vector quantities 
belong to the simultaneity submanifold. The main point, 
though, is that f a should in no sense depend on quan- 
tities associated with the future (or past) particle states. 
Taking p and U a together as constituting the local state 
of the particle, then the determination of f a at p should 
not depend on future or past values for these quantities, 
neither explicitly nor implicitly. 

Technically speaking, it is the F a p tensor, rather than 
A a itself, that is used in the construction of f a . As it 
happens, the F a p tensor does depend on future particle 
states, because this quantity is obtained from all partial 
derivatives of A a . Thus, it is not only A a at the point p 
itself that is consulted, but also, the value of A a at nearby 
points, displaced from p in all four spacetime directions. 
This situation might appear to violate Postulate 2, but in 
fact, it does not. The reason is that the full F a p tensor is 
not actually required to construct /" — only that portion 
of the tensor that is projected onto the local simultaneity 
submanifold is used. Consequently, although A a values 
at points other than p are required, these points occur 
simultaneously with p, lying neither in its past nor its fu- 
ture, and so Postulate 2 is indeed satisfied. Similar argu- 
ments will prove relevant, also, for quantum force fields, 
although these arise in an altogether different manner. 

A more precise statement of Postulate 2 also requires 
a clearer specification of the dynamical "evolution coor- 
dinate," denoted £. Under arbitrary reparametrizations, 
£ — > £'(£)> the forms of Eqs. (37) and (39) are not pre- 
served, and — depending on interpretation — "time" com- 
ponents of the resultant force vectors do in fact arise. 
However, these should be properly regarded as fictitious. 
The fact is that the metric tensor, together with the lo- 
cal simultaneity submanifolds, give rise to a natural local 
choice for £, in terms of which Postulate 2 may be ex- 
pected to hold. This choice is such that £ = const within 
the simultaneity submanifold, and d^/ds = const for dis- 
placements perpendicular to the simultaneity submani- 



fold (where ds 2 is the line element, defined in Sec. II A). 
This specification of the evolution coordinate £ may be 
applied to both the non-relativistic and relativistic cases 
considered above — giving rise to £ = t and £ = r, respec- 
tively (more generally, £ may also include affine trans- 
formations of the above forms). With these choices for 
£, we have already seen that Postulate 2 holds for both 
non-relativistic and relativistic classical mechanics. 

B. Arbitrary time reparametrizations, as applied 
to classical electromagnetism 

As the canonical example of a relativistic classical force 
field with the above velocity-dependent form, we consider 
the case of electromagnetism in some detail. This will 
turn out to be of particular benefit for the development 
of relativistic quantum force fields, in at least two ways. 
First, it provides an understanding of how potential en- 
ergy contributions, in general, should enter into the rel- 
ativistic Lagrangian. The second benefit emerges from a 
reworking of the electromagnetic example in terms of an 
arbitrary time reparametrization, i.e. r — > A(r), which is 
not usually seen in standard treatments. In the electro- 
magnetic context, this introduces unnecessary complex- 
ity, but it provides important insight into how to handle 
the relativistic quantum case — for which, as we have seen 
in Sec. Ill, ensemble time A and proper time r are incom- 
patible. 

Let us start with the simplest case of a free particle. 
Let x a (X) denote an arbitrary timelike path, described 
by the arbitrary parameter A. For fixed endpoints, the 
solution trajectory is that which maximizes the elapsed 
proper time, At — i.e., it is a geodesic. From Eq. (2), the 
elapsed proper time along any path is given by 

ff dT \ „ f l~ l dxa dxP ~ 

AT -j{dx) dX -jy-^^^x^x dX - (41) 

One standard choice for A is the proper time, t, itself. 
Making this choice, and multiplying Eq. (41) by the neg- 
ative rest energy, — mc 2 , the maximization of At is con- 
verted into the equivalent minimization of the action, 

S = J L[x°] dr = J (-mc^^j-^^x-xUr, (42) 

where x% denotes dx a /dr — U a , and L[x"] is the La- 
grangian for the parameter choice, A = r. 

Regardless of the particular choice of A, the form on 
the right hand side of Eq. (41) should be used in the 
Lagrangian. This is because it exhibits explicit depen- 
dences on x\ = dx a /d\ = V a , and thus gives rise to an 
Euler-Lagrange ordinary differential equation (ODE) for 
the solution trajectory, x a (X), in terms of the first and 
second A derivatives. Specifically, action minimization 
leads to the following Euler-Lagrange form: 

dL_ _ d_ ( dL\ 

dx^ d\\dx^J ~° ( 3) 
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For the standard choice A = r, and the flat Minkowski 
spacetime metric presumed here, the result is Eq. (39) 
(with f a =0). For curved spacetime manifolds, or flat 
manifolds described via curvilinear coordinates, A^ 1 , tak- 
ing A = r leads to the standard geodesic equation, 33 ' 43 

X» T + Y^X»X* = 0. (44) 

To a large extent, the above procedure is straightfor- 
ward because of the presumed A = r form, and the con- 
straint that that induces. Even in this context, there 
are some ambiguities that can arise. Specifically, the 
square root quantity in Eq. (42) is necessarily equal to 
one. Why not simply replace this expression with 1, 
or otherwise multiply or divide L[x%] by as many such 
square root factors as we wish? Likewise, incorporating 
the A = t constraint explicitly into the optimization us- 
ing the theory of Lagrange multipliers, one should be able 
to add or subtract additional terms of this form. Even 
for the free particle case, this can lead to trouble — e.g., 
to S — f(—mc 2 )dT, which is useless from the point of 
view of generating any ODE for x a (r), much less the cor- 
rect one. The situation is even more delicate, however, 
when a potential energy contribution is introduced into 
the Lagrangian. 

Consider, then, the standard form of the Lagrangian 
(in r) for a single relativistic particle with charge q, acted 
on by the electromagnetic vector potential, A a : 



L[x a ,x«] = (~ m c 2 )^-^ VaP dx? dx P T + ^A a (45) 

The potential energy contribution [i.e., the second term 
on the right hand side of Eq. (45)] introduces an explicit 
Lagrangian dependence on the coordinates x a , through 
the A a vector field. The Lagrangian form of Eq. (45) 
leads via Eq. (43) to the correct electromagnetic ODE 
i.e., Eq. (39), with 

f a = - c v ap F Pl xl and F M = dp Ay - d 1 A p . (46) 

The above derivation requires explicit substitution of 
the Eq. (2) equality — but only after the Lagrangian par- 
tial derivatives of Eq. (43) have been applied. Note that 
the one factor of the square root quantity in the first 
term of Eq. (45) is necessary to generate the kinetic en- 
ergy contribution in the resultant equation of motion. 
The potential energy term of the Lagrangian does not 
include this square root factor; if it did, the resultant 
ODE would be incorrect. Likewise, any additional fac- 
tors of the square root quantity in the first (kinetic) term 
of L would also lead to incorrect ODEs. 

We will apply these findings to relativistic quantum 
force fields, but first we must also consider the general 
A case explicitly. This can be regarded in terms of the 
coordinate transformation, r — > A(r). Since the action, 
S, must be a scalar invariant, L must transform as a 
scalar density of weight W = — l — leading at once to the 



following form for the Lagrangian in A: 



L[x a ,xf] = {-mc 2 )^-^^ x a x x^ + ^ x a x A a (47) 

At this point, we must distinguish between two possi- 
bilities: the case where the explicit relation A(r) is known 
a priori, vs. the case where it is not. In the former 
case, A(r) implies an explicit constraint, analogous to 
(but more complicated than) that of Eq. (2) for A = r. 
One can therefore apply a procedure similar to that used 
for A = r i.e., the explicit constraint is invoked only 
after the Lagrangian partial derivatives have been eval- 
uated. The resultant ODE is the general-A version of 
Eq. (39), i.e., 

mxl^rtr+m^xl, (48) 

where t\ = dr/dX, and f a is from Eq. (46). That 
this form is correct may be readily verified — e.g., for the 
choice A = t = x°/c, in terms of which the classical elec- 
tromagnetic equations of motion are very well known. 

Let us now imagine — as is the case in the quantum 
context — that the relation A(r) [or the inverse, r(A)] is 
not necessarily known a priori. This situation is evidently 
problematic, from the perspective of the A evolution of 
Eq. (48). At the initial time, it is straightforward enough 
to specify initial values for t\ and t\\ — or equivalcntly, 
A T and A TT . However, there are many functions A(r) that 
share the same initial conditions, any one of which might 
in principle be legitimately used; so how does the ODE 
"know" which choice to make, during the course of the 
propagation, if this is not already specified a priori? 

To address this question, the only recourse is to use the 
generic relation for t\, as implied by Eq. (41). Substi- 
tuting this form into Eq. (48) — or equivalently, deriving 
the Euler-Lagrange ODE from Eq. (47) without applying 
any explicit constraint — one obtains 

™x% x = (-^VPi4*l)r, (49) 

which is equivalent to Eq. (48) , but without the last term. 
Thus, t\\ = X TT = 0, or A(t) = At + B. In other words, 
this procedure automatically picks out an affine relation 
for A(t), if no relation is specified a priori. (The values 
of the constants A and B are uniquely determined from 
the initial values for A and A T , which are arbitrary). 

All of the above is perfectly consistent — albeit unneces- 
sarily complicated — in the context of relativistic electro- 
dynamics. The situation does not bode so well, however, 
for the relativistic quantum dynamical case. Here, we 
know that A and r do not satisfy any global functional 
relationship of the form A = A(r) — much less an affine 
relation. Locally — and along a single given trajectory 
of the ensemble — such a relation can of course be estab- 
lished, although it is not necessarily known a priori, even 
in this context. All of this suggests that, in the quantum 
case, we must proceed with caution. 
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Non-relativistic quantum dynamical 
considerations 



The trajectory-based formulation for non-relativistic 
quantum dynamics is, by this point, fairly well 
established. 20-26 Here, we review only those features that 
are most relevant for the relativistic generalization of 
Sec. V. Also, whereas our previous work has expressed 
the requisite quantities in terms of the Jacobian matrices, 
j\ and K\, here, we adopt a tensorial, generally covari- 
ant viewpoint, relying on the (spatial) metric tensor 7^, 
and its determinant, 7. 

The non-relativistic formulation makes use of a spatial 
coordinate transformation from x l to C 4 , for which the 
time coordinate, t, is unaffected. As in the relativistic 
case, the C % serve as trajectory labeling coordinates for a 
trajectory (or path) ensemble. As discussed in Sec. IV A, 
the simultaneity submanifolds are defined by surfaces of 
constant t in the 4d non-relativistic spacetime manifold, 
M. The x l — > C % transformation is thus parametrized by 
t, with a complete solution trajectory ensemble taking 
the form x l (t, C). In the x l coordinate frame, the metric 
tensor is presumed to be the usual Euclidean one, 5i m , 
whereas in the C l coordinates, the metric tensor is de- 
noted jij. The latter is a function of both space (C) and 
time (t). 

As in the relativistic case, Postulate 1 is presumed (or, 
it can be regarded as a logical consequence of Bohmian 
mechanics 7-13 ). The scalar probability density in C 
space, /(C), is thus independent of t, and Eq. (19) is 
presumed to hold. From Eq. (18), we have 



R(t,x) 



/(C) 1 / 2 



yl/4 



(50) 



where R(t, x) = p(t,x.y^ 2 is the usual wavefunction am- 
plitude, |^| — i.e., the square root of the usual spatial 
probability density, p(t,x) = The latter is a 3d 

density akin to j°/c from Sec. HID, but in the non- 
relativistic context, it can also be interpreted as the 4d 
scalar probability density of Sec. IIIC. 

We expect a scalar invariant action, S, of the form 



S = 



JjC[x l ,x l ,diX l ,...]d s Cdt 
= I L[x l ,x l ,diX l ,...]f{C)d 3 Cdt, (51) 



where x l = (dx l /dt)\C and di = d/dC % . The quantity 
C is the Lagrangian density (a spatial scalar density of 
weight W — —1), whereas L is a scalar invariant quantity 
(with units of energy) that we refer to simply as "the 
Lagrangian" (although technically, that term should be 
applied to the d 3 C integral of C.) The Lagrangian is a 
true scalar invariant (weight W = 0). In non-relativistic 
quantum mechanics, it takes the form 



L = 



-mS lm x l x m ~V{s) 



(52) 



8m 



/(C) 



7 



1/2 
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/(C) 



7 



1/2 



d. 



/(C) 



7 



1/2 



The metric tensor 7 depends explicitly on the <9jx'; the 
specific Lagrangian of Eq. (52) is therefore second order 
in the C l derivatives of x l (although an alternate, higher 
order choice of L will also be considered). 

The last term in Eq. (52) above is the (scalar invari- 
ant) quantum contribution to the Lagrangian, denoted 
— Lq. This term accounts for all interaction or "commu- 
nication" across trajectories — and hence, for all quantum 
dynamical effects. 21 ' 23,25 It can be expressed in a variety 
of ways, with the Eq. (52) form above being particularly 
explicit. Note that Eqs. (22), (32), and (50) imply that 
the quantity in square brackets can be interpreted as the 
spatial scalar invariant probability density, /* — or equiv- 
alently, as the usual p(t,x) = i?(i,x) 2 . Because /* is a 
scalar invariant, the partial derivatives di in Eq. (52) may 
be replaced with the corresponding covariant derivatives, 
Vi. Also, since (V*/)// is a true covariant vector (weight 
W = 0), f* may be replaced with /, if desired. Thus are 
we led to various alternate forms for Lq, e.g.: 



J Q 



8m 



Mi,x)rV^,x)]d>(t,x)] 



n 2 [vv(c)] [V,/(C)] 



8m 



/(C) 2 



(53) 



The expressions above are all manifestly covariant with 
respect to arbitrary coordinate transformations of the 
C % that do not depend (even parametrically) on t — i.e., 
transformations of the form C l — > C l = C l (C). Thus, 
for example, the uniformizing choice C = P gives rise to 
these same expressions, but with C replaced with P, and 
/(C) replaced with 1. It is useful, however, to extend the 
range of coordinate transformations to include those that 
do depend on i, at least in the parametric sense indicated 
above. We refer to such general coordinates as X 7 ', and 
see that they include both C l and x l as special cases. In 
Xi coordinates, the generalized scalar probability den- 
sity (weight W = -l) is denoted p(t,X) = R(t,X) 2 . We 
can therefore write the general form of Lq as 



Lc 



h 2 [VV(t,X)] [V jP (t,X)} 



8m 



p(t,X) 2 



h 2 [V j R(t,X)] [VjR(t,X) 
2m R(t,X) 2 



(54) 



Consider the second form in Eq. (54) above, for the 
specific choice A- 7 = x . The covariant derivatives 
become ordinary partial derivatives (gradients), even 
though R is a scalar density of nonzero weight (W = 
-1/2). The resultant 



C Q = (n 2 /2m)[Vi?(i,x)]-[Vi?(t,x)] 



(55) 



closely resembles the gradient contribution to the KG 
Lagrangian — but with the partial derivative components 
restricted to the simultaneity submanifolds (Sec. VE). 



15 



Another important dynamical quantity is the quantum 
potential, Q, which, in generally covariant form, is given 
as 



Q 



h 2 V J Vji?(i,X) 
~2m R(t,X) 



(56) 



The quantum potential is a true scalar invariant quantity. 
For A J = x l , V J Vj becomes the usual Laplacian, and so 
Eq. (56) reduces to the familiar expression from Bohmian 
mechanics. 8-11 In this form, 2RQ is the contribution to 
the Euler- Lagrange PDE that results from £q, treating 
R (rather than x) as the dependent field. This matches 
the corresponding KG PDE contribution — but again — 
restricted to the simultaneity submanifolds (Sec. VE). 

There is another interesting relation between Q and 
Lq, that arises in a — C % coordinate representa- 
tion. Here, the quantity R(t 7 X) in Eq. (56) becomes 
/(C) 1 / 2 . However, it is more convenient to replace / 
with the scalar invariant /*, as discussed above. Finally, 
since the covariant Laplacian V J Vj = VjV- 7 is now be- 
ing applied to a true scalar invariant, it can be expressed 
explicitly in terms of ordinary partial derivatives using 
the Laplace-Bcltrami form. The result is: 



Ql 



X" C i,X C i C j , x l C i C j Ck ] 



1 



di 



2m V7 1/4 / 1/2 



(57) 



' \ 7 l/4 



where x ci = diX 1 — dx l /dC l , etc. 

Equation (57) is the trajectory-based form of Q, which 
is appropriate when solving for the solution trajectory 
ensemble, x (t, C). Since 7 depends explicitly on x ct , as 
discussed, the quantity Q is third order in the C n deriva- 
tives of x l . The relation to Lq is that one may substitute 
the Lq form used in the Lagrangian of Eq. (52), with 
Lq = Q, without altering the resultant Euler-Lagrange 
PDE for x l (t, C). This change therefore amounts to a 
change of gauge — with the Lq = Q choice offering cer- 
tain advantages, despite the higher order, as discussed in 
previous work. 22 ' 24 ' 25 In any case, it can be shown that 
Eq. (57) is equivalent to the J l i and K\-based expressions 
for Q that have been derived previously. 

Whether or not Q is used in the Lagrangian, it invari- 
ably plays a direct role in the resultant Euler-Lagrange 
PDE for x l (t, C). Specifically, the covariant quantum 
force vector is the (negative) gradient of the quantum 
potential — i.e., 



Q 



(58) 



in generally covariant form. The quantum force enters 
the PDE in the manner expected: 



mx l + 5 



lmdV(x) 



dx r 



f l Q 



2m 



where 



6 Lm K" m d k 



7 V2yj q 



/V 2> 



M y/4_ 



and K k m = dC k /dx m . 

The Euler-Lagrange PDE for x l (t, C)— i.e., Eq. (59)— 
is second order in time (t) and fourth order in space (C). 
Note, however, that no time derivatives enter into the 
expressions for the quantum (and classical) forces. At 
a given point p, therefore, the /q vector belongs to the 
orthogonal subspace of T p , and is otherwise only con- 
structed from quantities that also belong to the simul- 
taneity submanifold associated with p. Postulate 2 is 
therefore satisfied, for the trajectory-based formulation 
of non-relativistic quantum mechanics. 



V. RELATIVISTIC QUANTUM DYNAMICS 

A. Relativistic quantum dynamical equation of 
motion 

We have taken considerable lengths to lay the ground- 
work for our main objective: an exact, self-contained, 
trajectory-ensemble-based PDE, describing the relativis- 
tic quantum dynamics of a spin-zero free particle. This 
foundational effort has been well spent, in that it gives 
rise to an essentially unique relativistic quantum formula- 
tion, requiring only the barest minimum of assumptions. 

Let us start with the relativistic quantum potential, 
Q. From previous non-relativistic work, it is known that 
the allowed trajectory-based form for Q is essentially 
unique. 22,24 ' 25 The same basic form should be required 
in the relativistic context — although here, there might in 
principle be some question as to which metric tensor (i.e., 
3d or 4d) and/or probability density quantity should be 
employed. Postulate 2, along with the broader discussion 
of Sec. IV A, completely alleviates this ambiguity. In par- 
ticular, since Q directly gives rise to quantum forces, it 
must be constructed on the relativistic simultaneity sub- 
manifold, using the 3d relativistic spatial metric tensor 
7 of Eq. (13) and Sec. Ill A. (More technically — and 
for completely general coordinates X^ and metric ten- 
sors g — the spatial metric tensor 7 that is used for this 
purpose is the induced metric tensor on the simultaneity 
submanifold, 33,43 i.e. the pullback of g). Likewise, the 
relevant probability quantity is the spatial scalar prob- 
ability density, /(C). As discussed in Sec. IV B, /(C) 
"lives" on the simultaneity submanifolds, and, is there- 
fore unaffected by reparametrizations, A — > A' (A), of the 
ensemble time coordinate, A. 

Without further ado, then, we posit a relativistic quan- 
tum potential of the form of Eq. (57), i.e. 



(59) Q[xc',x^, iC] ,x^ tCjck } - ^ m 



(60) 



r 1 / V di 



7 7 • "i \ .,1/4 
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One important difference between Eqs. (57) and (60) is 
that in the latter, t = x° jc is a coordinate, rather than 
a parameter. Thus, for example, dx° /dC l appears ex- 
plicitly in the expression for 7, in the relativistic case 
(along with dx l jdC % , of course). Another important dif- 
ference is that the di = d/dC l partial derivatives occur 
at constant A, rather than at constant t. We might have 
expected the relativistic partial derivatives to be at con- 
stant r; the fact that these occur at constant A, instead, 
is very significant. In any event, the relativistic quantum 
potential as defined by Eq. (60) is invariant with respect 
to reparametrizations of both A and C. 

Another significant feature of Eq. (60) is that the ex- 
plicit partial derivatives, di, extend over only the spa- 
tial, C l components of the natural coordinates [X^ from 
Eq. (12)]. The lack of a timelike, A contribution in effect 
implies that the 4d Laplace-Bcltrami operator has been 
projected onto the simultaneity submanifold — a neces- 
sary consequence of Postulate 2, of course. It is possible 
to develop a completely general covariant expression for 
Q, in terms of completely arbitrary spacetime coordi- 
nates. In this case, however, the relativistic analog of 
Eq. (56) must be substantially modified, to explicitly in- 
corporate the requisite projection tensors associated with 
the pullback of g. 33 > 43 We do not find it profitable to do 
so here. 

The relativistic quantum force is obtained as the (neg- 
ative) covariant gradient of the quantum potential scalar 
invariant of Eq. (60). Again, Postulate 2 requires that 
only the gradient components corresponding to the C % 
coordinates be considered — i.e., the gradient is restricted 
to the simultaneity submanifold. In a completely gen- 
eral coordinate system, each of the gradient components 
of the force vector must be explicitly projected onto the 
submanifold, thus modifying the form of Eq. (58). In the 
x a inertial coordinates, for example, this projection is ef- 
fected using the twelve \x = i components of the inverse 
Jacobi matrix — i.e., the K' l a — dC l /dx a . The inertial 
components of the relativistic quantum force vector, fg, 
are thus given explicitly by 



fc 



^K k fj d k Q 
2m 



(61) 
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As discussed in Sec. IV A, the time coordinate associ- 
ated with the dynamical evolution is r, rather than A. 
For the relativistic quantum free particle, the only dy- 
namical force is the quantum force, defined by Eq. (61). 
We therefore expect the relativistic quantum equation of 
motion to be given by Eq. (39), with f a = fg- This equa- 
tion thus constitutes a new (to our knowledge) dynam- 
ical prediction, that can in principle be tested against 
experiment. It can be interpreted as a PDE, but one 
that is not equivalent to the KG PDE (Sec. VE). The 
proposed dynamical law is unambiguous, explicit, and 
exact. However, in its present form, it has a "mixed" 
character, whereby A is the time coordinate used in the 



determination of Q and /q, but r is what is used in the 
dynamical law. Later, we will derive a single, consistent 
form, expressed solely in terms of A. First, however, we 
will demonstrate that the proposed dynamical law is in 
fact the Euler-Lagrange equation that results from ex- 
tremizing the corresponding relativistic quantum action. 



B. Extremization of the relativistic quantum action 

The relativistic quantum dynamical law proposed in 
Sec. VA [i.e. Eqs. (39) and (61)] should be derivable as 
the Euler-Lagrange PDE obtained from the extremiza- 
tion of some suitable action quantity. Moreover, from a 
perspective of general covariance, this action should be a 
true scalar invariant, obtained as an integral over d 4 X, 
for arbitrary coordinates, X^. In analogy with Eq. (51), 
therefore, we expect 



S = I c[x a ,x% .... d'\\ 
L[x a , x'xn , . . . 



(62) 



where p(X^) is the scalar probability density of Eq. (24), 
and the Lagrangian, L, is a true scalar invariant with 
units of energy. 

If X^ is taken to be a set of natural coordinates, then 
Eq. (25) implies that Eq. (62) receives an extra factor 
°f V~ .9oo = t\ in comparison with Eq. (51), where t\ 
is defined as in Sec. IV B. This fundamental difference 
from the non-relativistic case is associated with the fact 
that it is t rather than A that is the dynamical evolution 
coordinate in the relativistic context. In any case, from 
the discussion in Sec. IV B [e.g., Eqs. (41) and (47)] and 
Sec. IV C [Eq. (52)] the explicit form of the relativistic 
quantum action, in natural coordinates, is as follows: 



- J (-mc 2 )^ 



(-mc 2 )\/-^T] a p x* x% - L Q 



T X f(C)d 3 Cd\ 



J (— mc 2 ) 
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f(C)d 3 Cd\ (63) 



Since t\ depends explicitly on the xjr^ and implicitly 
on the X^, its presence in the Lagrangian density defi- 
nitely complicates matters, vis-a-vis the resultant Euler- 
Lagrange PDE. Matters are further complicated by the 
fact that the A coordinate is not defined a priori, as dis- 
cussed in Sec. IV B. We therefore anticipate difficulties, 
or at least complications, in any Euler-Lagrange deriva- 
tion based on natural coordinates, X^ = (cA, C). The 
most expedient way around this difficulty is to adopt a 
kind of hybrid approach, involving a change of coordi- 
nates, X» ->• Y U {X^) = (ct(V^),C), where the time 
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coordinate for the new frame, Y u , has been changed from 
A to r. 

An advantage of working in the Y v frame is that the 
time coordinate r is now well determined. Overall, how- 
ever, the above change of coordinates does not necessarily 
constitute an improvement; for example, the new metric 
tensor, g, is no longer block diagonal. Also, the restric- 
tion to the simultaneity submanifolds of quantities such 
as Lq and Q would require explicit use of projection ten- 
sors. On the other hand, we recall Eq. (9), which allows 
us to redefine the proper time coordinate, r, by applying 
an arbitrary, path-dependent shift, A(C). This enables 
us to fix the particular choice of t, such that one of its 
contours, at least, coincides exactly with a simultaneity 
submanifold. 

Let us, then, hereby choose A(C) such that A = cor- 
responds to r = 0. As a consequence, for all points p 
belonging to the specific simultaneity submanifold speci- 
fied by A = r = 0, the partial derivatives c*j|A and di\r are 
the same, and so the latter can replace the former in the 
expressions for Lq, Q, and fg. Along this simultaneity 
submanifold, therefore, the Lagrangian in Y v becomes 



(-mc 2 ) 

l / h 

8 \ mc 
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dx a dxf 3 
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where the Y % = C l partial derivatives, di — explicit in 
Eq. (64) above, and implicit in the definition of the Xy v 
quantities used to construct the spatial metric tensor 7 — 
are taken at constant r, rather than at constant A. 

The advantage of Eq. (64) is that it is now very 
straightforward to derive the corresponding Eulcr- 
Lagrangc PDE — at least as this is restricted to the A = 
t = simultaneity submanifold. The result is indeed the 
dynamical law of Sec. V A — i.e., Eqs. (39) and (61), with 
f a = Jq. Note also that these equations have no ex- 
plicit dependence on the r coordinate itself — only on the 
differential, dr, which is invariant with respect to trans- 
formations of the form of Eq. (9). This is significant, 
because it implies that the dynamical law must hold not 
only for the A = r = simultaneity submanifold, but for 
all simultaneity submanifolds — and therefore, across the 
entire spacetime manifold, M. 



C. Conversion to natural coordinates, and 
"dynamical" ensemble proper time, T 

Our present theory encompasses a hybrid dynamical 
equation, wherein r is the time evolution coordinate as- 
sociated with the quantum force vector f a = /q, but 
the latter is itself obtained at constant values of the en- 
semble time coordinate, A. (We henceforth drop the Q 
subscript on the quantum force vector, f a ). In practi- 
cal terms, such an equation is not so useful for actually 



effecting propagation of the solution trajectory ensem- 
ble, x a (X^). We therefore seek to convert the dynamical 
equation to a more consistent form involving only the 
natural coordinates, X^ = (cA, C), that can be directly 
integrated with respect to A, and without recourse to t. 

Since the equation of motion is second order in time, 
the quantities to be propagated may be taken to be 
x a (C) and U a (C). We seek new time evolution equa- 
tions for the derivatives of these quantities with respect 
to A. From Eqs. (25) and (27), and also Eq. (39), these 
are found to be: 
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~d~X 


~~ ° dX° 
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dU a 


dU a 
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= r x U c 
f° 



m 



(65) 



We therefore have what we need in Eq. (65) — provided 
that an explicit expression for t\ can be obtained. In 
general, this quantity varies across spacetime, and even 
across a simultaneity submanifold. In addition, it will 
change under a reparametrization of the A — > A' (A) form. 
Also, in the special case of quantum inertial motion — 
provisionally defined via f a — for all x a (a more pre- 
cise definition will soon be introduced) — an arbitrary A 
should satisfy Tx(X^) = h(X), for some function h(X). 
In other words, for the inertial case, t\ should not vary 
across a given simultaneity submanifold. Choosing A to 
be an ensemble proper time, T, one has h(T) = 1 for 
the inertial case, by definition. In any case, constant t\ 
across a given simultaneity submanifold implies f a = 0, 
and so conversely, we expect nonzero quantum forces to 
give rise to local variations in t\. 

What precise relationship should one expect to see be- 
tween t\ variations and the quantum force? Consider 
that a A reparametrization has the effect of rescaling t\ 
uniformly across a given simultaneity submanifold. As 
discussed, such a reparametrization has no effect on f a . 
Therefore, regardless of how A defined, it is the relative 
changes in t\ — i.e., di(T\)/r\ — that should appear in the 
quantum force relation. 

In any event, an explicit form for this relation may be 
obtained upon transforming Eq. (65) to a natural coor- 
dinate frame. Treating t in this context as a parameter, 
or scalar field, we have already seen that both x a and 
U a are true four-vectors. The covariant derivative of the 
U a vector field, i.e. 



(66) 



is also a true tensor invariant, where the second equation 
above is the generally covariant form. Acting on an ar- 
bitrary displacement vector, dx^ , this tensor yields the 
corresponding change in the velocity four- vector, dll a . 
To compute the change in U a along the direction of mo- 
tion itself, Tp a should be contracted with U@/c. Upon 
dividing both sides by t\, the second part of Eq. (65) 
thus becomes 



^ = ueT fi ° = r 

dr m 



(67) 
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WTS = u v v v u^ = u v (djji* + r^j/' T ) 



where the second equation above is the generally covari- 
ant form. 

From Eqs. (4) and (20), the natural coordinate form 
of Z7 M is easily found to be 



UP- 



,0,0,0 



(68) 



From Eq. (68), it is also easy to verify the natural co- 
ordinate generalization of the normalization condition, 
Eq. (3). Likewise, the orthogonality of U a and f& — i.e., 
Eq. (40) — must also hold in the natural coordinate frame: 



(69) 



From Eqs. (68) and (69), and the block-diagonal nature 
of g [Eq. (13)], it follows that 



(70) 



in a natural coordinate representation. The fact that 
f° = is consistent with the requirement that / M live 
on the simultaneity submanifold. However, the metric 
ensures that even though the / M vector is always "hor- 
izontal," the vector is nevertheless always "vertical" 
(Fig. 3). 

The results of the previous paragraph must be recon- 
ciled with Eq. (67). In natural coordinates, 

r S „ = I s »»a„ = i(^)M (71) 

Substituting Eqs. (68) and (71) into Eq. (67) thus yields 
/ = 0, as expected. It is the expression for the spatial 
components of force, however, that yields the desired t\ 
relation. Here, the result 



r i 00 = -li ij d j g 00 = ^d j T 2 x 



gives rise to 



m \ t x 



(72) 



(73) 



which is exactly of the form desired, as discussed above. 

Equation (73) above may be rewritten in more general 
and suggestive fashion as 



,f = mc 2 ^ d 3 [\og(A(X)T X )} , 



(74) 



where A(X) is constant across a given simultaneity sub- 
manifold, but depends arbitrarily on A. Clearly, the func- 
tion A(X) defines a particular A parametrization; it has 
units of A over time. The most natural choice is simply 
A(X) = 1. With this choice, A has units of time, and sat- 
isfies T\ = 1 in the case of inertial motion. It is therefore 
an ensemble proper time, A = T. We call this choice the 



"dynamical ensemble proper time"; henceforth, it is the 
only choice for T that will be considered. 

It is necessary to demonstrate consistency of Eq. (74) 
across all values of i. This follows easily from the fact 
that /' is the (spatial) gradient of the quantum potential 
[Eq. (58)], thus enabling a direct comparison between t\ 
and Q. Choosing A = T, this leads to the fundamental 
result, 



dr_ 
df 



Tf = exp 



Q 



mc* 



(75) 



The importance of Eq. (75) cannot be overstated. In 
addition to providing exactly the t\ — tj- form needed 
to propagate Eq. (65), this expression is imbued with 
fundamental physical significance. Specifically, it states 
that the quantum potential induces a time dilation or 
contraction for individual trajectories, much like gravi- 
tational time dilation. Also like the case of gravity — but 
unlike most other physical forces — the absolute quan- 
tum potential matters. Thus, under the transformation, 
Q — > Q + const, there is a corresponding change in T— 
a constant rescaling, in fact — even though the quantum 
forces and accelerations are unaltered. Put another way, 
Eq. (75) requires that quantum inertial motion satisfy a 
more restrictive condition than stated above; specifically, 
we must have Q = everywhere, rather than just /'=0. 

An important difference from the case of gravity is that 
the quantum potential Q may be either positive or neg- 
ative. In other words, the passage of proper time for a 
given quantum trajectory may be slower or faster than 
that for the whole ensemble, depending on the sign of Q. 
There is also a curious inverse behavior, whereby lower 
(i.e. more negative) values of the quantum potential lead 
to faster evolution of t (relative to T) exactly the op- 
posite of the gravitational case. More specifically, in the 
weak field limit, one has 
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.9oo 
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(76) 



where the left hand side of Eq. (76) is for GR, and the 
right hand side is for the present quantum theory. Thus, 
it is —Q, rather than Q itself, that plays the role of the 
gravitational potential energy, m$. 

On the other hand, the above behavior makes perfect 
sense, from another perspective. In standard Bohmian 
mechanics, regions of space where Q is positive may be 
(loosely) identified as classically allowed regions, whereas 
negative Q generally denotes classically forbidden regions 
into which quantum tunneling can occur. 10,11 In the clas- 
sically allowed regions, therefore, the proper time for a 
given trajectory is dilated — as would be predicted by 
classical relativity theory. Quantum tunneling, on the 
other hand, speeds up the passage of proper time — a phe- 
nomenon with no classical relativistic analog. Although 
the present work focuses on the flat spacetime manifolds 
of SR, the features of relativistic quantum dynamics dis- 
cussed above suggest an extremely interesting interplay 
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between gravitational and quantum forces, in the context 
of quantum gravity. 

In any event, substitution of Eq. (75) into Eq. (65) pro- 
vides us with our goal for this subsection — a set of dy- 
namical equations, consistently expressed in terms of any 
choice of natural coordinates, = (c A, C), that can be 
directly integrated over A. [To obtain the general-A form 
from the A = T form, one simply introduces a factor of 
A(X) into the exponent of Eq. (75).] These equations, 
for dx a /d\ and dU a /d\, are first order in A and fourth 
order in C, with all quantities on the right hand sides 
of Eq. (65) obtained from a single simultaneity submani- 
fold [i.e., from /(C), x Q (C), and various di partial deriva- 
tives, all at fixed A] . We also find it convenient to express 
this dynamical law as a second order (in time) PDE, using 
the specific choice of natural coordinates, = (cT, P), 
where T is the dynamical ensemble proper time, and P 
are the uniformizing coordinates of Sec. HIE. The result 
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In practice, the last form of the last line above may 
be more convenient than the middle form, because it in- 
volves explicit partial derivatives of the x a rather than 
of the X^ (i.e., the elements of the Jacobian matrix J a ^ 1 
rather than the inverse Jacobian K^ a ). It can be easily 
derived, simply by raising the index of the force vector 
prior to transforming from natural to inertial coordinates 
[e.g., via Eq. (73)], and exploiting Eqs. (13) and (70). It 
can also be shown that the various constraints imposed 
on U a , f a , etc., as discussed previously in this document, 
are all preserved under the evolution of Eq. (77). 



D. Classical, inertial, and non-relativistic limits 

As a necessary check on the viability of our candi- 
date relativistic quantum dynamical law, as presented in 
Sees. VA and VC, we must consider various asymptotic 
limits for which the correct behavior is known. Two ob- 
vious limits that must be considered are the relativistic 
classical limit, and the non-relativistic quantum limit. 
Fortunately, the trajectory-based formulation adopted 
here renders both of these comparisons very straightfor- 
ward. In addition to these two limits, one should also 
consider the limit of quantum inertial motion, alluded to 
several times previously in this document. 

The relativistic classical limit is the most straightfor- 
ward. Consider the full, relativistic quantum action, as 
presented in Eq. (63). We define the "relativistic classi- 
cal limit" of this action as what ensues by setting % = 0. 
Of course, the second, Lq term vanishes completely in 



this case. One is indeed left with the relativistic clas- 
sical action of Eq. (42), except that it is integrated by 
/(C) d 3 C. This integration is immaterial, however, be- 
cause the relativistic classical Lagrangian density — being 
truly classical — does not depend explicitly or implicitly 
on the x^i. Thus, in the relativistic classical limit, 
each trajectory of the solution ensemble propagates com- 
pletely independently of the others, and in accord with 
the relativistic classical equation of motion, Eq. (39). 
This behavior is completely as expected — being analo- 
gous, e.g., to the Hamilton- Jacobi trajectory ensemble in 
the non-relativistic classical theory. 

Note that in both relativistic and non-relativistic con- 
texts, individual trajectories may cross each other, in the 
classical limit. For semiclassical and related approxima- 
tions to non-relativistic quantum mechanics that operate 
in the classical limit, trajectory crossings lead to well- 
known headaches such as caustics, divergent probability 
densities, multivalued field functions, etc. For such ap- 
proximations, also, it is well-known that "as h goes to 
zero" is a rather imprecise phrase, whose specific mathe- 
matical interpretation must be handled very delicately. 
In the trajectory-based formulation of non-relativistic 
quantum mechanics, for example, setting Ti — leads to 
the semiclassical or Hamilton- Jacobi behavior described 
above (e.g., with crossing trajectories) — whereas setting 
% to any finite value, no matter how small, leads to a 
qualitatively very different ensemble with no trajectory 
crossings. 10 ' 11 

The same is true for the present relativistic quantum 
formulation. Thus, it is necessary to distinguish the case 
where h = identically (considered above) from that 
where % approaches zero (in a certain well-defined sense 
given below). Simply put, in the latter case, there is 
still correlation across the trajectories in the ensemble, 
even though quantum dynamical effects become arbitrar- 
ily small. It is easy to show that this limit corresponds 
to the special case of quantum inertial motion, for which 
X^ = (cT, P') -> a;"'— with P' a rescaled version of P, 
and x a some inertial coordinate frame [sec Eq. (5)]. 

More precisely, quantum inertial motion is defined via 
the limit, 



mc 2 
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mc 
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djpjt, x) 
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(78) 



The physical significance of Eq. (78) is clear; it says that 
quantum inertial behavior emerges when the length scale 
of the probability distribution is extremely large com- 
pared to the Compton wavelength, h/mc. Thus, p(t,x) 
approaches "infinite broadness," or effectively uniform 
behavior. There is an obvious analogy with the Ein- 
stein Equivalence Principle, except that in the relativis- 
tic quantum context, "sufficiently small regions of space- 
time" must still be larger than the Compton wavelength. 

One can equally well define the quantum inertial limit 
via (Q/mc 2 ) — > 0, from which Eq. (75) leads to tj — > 1 or 
T — > t — a necessary condition for quantum inertial mo- 
tion. Because the quantum potential approaches zero, 
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the same must be true of the quantum force vector, /q, 
implying that the quantum trajectories follow nearly per- 
fect straight-line orbits, as in the classical limit case. 
However, unlike the classical limit case, the trajectories 
within a given ensemble are all parallel, in accord with 
Eq. (7). This assertion is easy to prove by contradic- 
tion. Suppose that two trajectories in the ensemble are 
not parallel; they must cross at some point p in M. This 
violates the presumed ensemble conditions of Sec. II, but 
according to Postulate 1, the probability density p must 
diverge at p, which also violates Eq. (78). Thus, in a 
global sense at least (i.e., across all M), quantum iner- 
tial motion implies Eq. (7). 

Finally, we consider the non-relativistic quantum limit. 
Here, we may expect that Eqs. (63) and (77) will reduce 
to the corresponding expressions for the non-relativistic 
trajectory-based formulation presented in Sec. IV C — but 
with V(x) = 0, as is appropriate for a free particle. First, 
we define the relativistic parameter (3 (not to be confused 
with the index (3) in the usual fashion, as the speed of 
the particle in units of c: 



x l Xl 

c c 



(79) 



where U a = (U°,U l ) and x l = dx l /dt. By definition, 
< (3 < 1, with ^40 taken to be the non-relativistic 
limit. 

To address this limit, it is therefore appropriate to ex- 
pand all quantities in Eq. (63) in powers of (3, in which 
context it is convenient to consider the classical (first 
term) and quantum (second term) contributions sepa- 
rately. The classical contribution has no explicit depen- 
dence on C, nor on its partial derivatives; we may there- 
fore take A to be any trajectory parametrization, with the 
choice \ — t a particularly convenient one. This yields 



J(-mc 2 )^/l-^f(C)d 3 Cdt 



(80) 



J (-mc 2 + ^mx l x^j /(C) d 3 C dt, 



where the right hand side above is the expansion of the 
left hand side, to second order in (3. The zeroth or- 
der contribution is the usual rest energy, which — being 
a constant — has no effect on the dynamics, and can be 
ignored. There is no first order contribution. The second 
order contribution is the expected non-relativistic classi- 
cal kinetic energy, i.e., the first term in Eq. (52). 

To maintain f3 order consistency, we should expand the 
quantum contribution to Eq. (63) to second order in (3 as 
well. Let At and Ax denote roughly the range of varia- 
tion of these respective quantities, for a given trajectory. 
In the relativistic limit of small /3, we have At s=y At, and 
so the quantum force / Q is on the order of mx a . This 
leads to 
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where "w" is to be interpreted here as meaning "on the 
order of." Similar arguments may be used to show that 
Q w Lq. Thus, the quantum contribution to Eq. (63) 
is (implicitly) second order in (3. We need therefore only 
consider the lowest order contribution in the explicit (3 
expansion of this quantity, which simplifies matters con- 
siderably. 

To zeroth order in (3, the trajectory is at rest, and so 
t w r, U° w c and U l w 0. The T and t contours are 
identical to zeroth order; likewise, the intersections of 
the C l contours (which define the trajectory) and the x l 
contours are identical. As a consequence, (d/dC l )\j- w 
(d/dC l )\ t - To lowest order, the J a M tensor thus becomes 
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From Eqs. (10) and (82), and the discussion in Sec. V A, 
the spatial metric tensor becomes 



^ij — J iJ^j ~ ^lr. 



dx 
dC* 



dx m 



(83) 



where the final expression above is the lowest order ap- 
proximation. But this is identical to the non-relativistic 
jij tensor of Eq. (51). Finally, from Eqs. (75) and (81), 
we have tt ~ 1 at lowest order. Thus, dT ~ dr rss dt, 
and we can replace dX in the integration of the second 
term in Eq. (63) with dt. Combining this with the second 
order result for the first term as derived above, one ob- 
tains exactly the non-relativistic Lagrangian of Eq. (51) 
[with V{x) = 0]. 

In principle, the above Lagrangian should lead, via the 
usual Euler-Lagrange procedure, to a set of dynamical 
PDEs that are consistent with the (3 — > limit of the rela- 
tivistic quantum equation of motion. Moreover, since the 
zeroth order (rest energy) contribution to the relativis- 
tic quantum Lagrangian is a constant, and since there is 
no first order contribution, we should need only consider 
the lowest order contribution to the relativistic quantum 
PDE, Eq. (77) (with a = I). To lowest order, the last 
term in this equation is zero, because dx 1 /dT w c Jo ~ 0. 
Likewise, the exponential factor is unity, to lowest order, 
and the partial derivatives with respect to T on the left 
hand side may be replaced with t derivatives. The re- 
sult, therefore, is indeed identical to Eq. (59) [again with 

y(x) = o]. 



E. Comparison with Klein-Gordon equation 

The standard description for a massive, spin-zero, rela- 
tivistic quantum particle is, of course, that of the famous 
Klein-Gordon (KG) equation. 3 ' 5 - 10,31_37 This is a wave 
PDE, for which the independent variables are x a , and 
the dependent quantity is the real- or complex-valued 
wave field, $(x Q ). The KG PDE takes the form 



p 2 <S> - d a d a $ = 0, 



(84) 
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where fi is a positive real constant (not to be confused 
with the index). The KG PDE is linear, and second or- 
der, in both the spatial coordinates, x , and the time 
coordinate, x°. It can be derived from a simple KG La- 
grangian density, involving only the first order derivative 
quantities, <9 Q $ = <& x a, as well as the wave field, $, itself: 



£[$,a Q $] = [p 2 <t> 2 + (d a <f>)(8 a $)] 



(85) 



The solution wave field — i.e., a $>(x a ) that satisfies 
Eq. (84) — extremizes the integral of Eq. (85) over d 4 x 
(or equivalently, over d 3 xdt). Because the KG PDE is 
linear, the units and magnitude of $ are immaterial. 

In non-relativistic quantum mechanics, one can estab- 
lish an equivalence between the trajectory-based nonlin- 
ear PDE of Eq. (59), describing the solution trajectory 
ensemble x l (t,C), and the linear TDSE PDE, describ- 
ing the complex- valued wavefunction field, ^(t, x). It is 
therefore natural to consider whether or not a similar 
equivalence can be established in the relativistic quan- 
tum case — perhaps leading to the KG PDE, or perhaps 
to a different wave PDE. To this end, we seek to con- 
vert the trajectory-based relativistic quantum action of 
Eq. (63) into a form involving only density fields, rather 
than explicit x\^ quantities. 

Using the constraint of Eq. (2), we first replace the 
square root quantity in the first line of Eq. (63) with a 
factor of unity. Such a modification is necessary in order 
to remove the explicit dependences; however, from 
the discussion in Sec. IV B, this change may lead to a fun- 
damentally different Euler-Lagrangc PDE. Next, Eq. (25) 
is used to replace r A /(C) with p(X») = R(X^) 2 , and 
Eqs. (22), (23), (24), and (54) are applied in order to 
rewrite Lq in terms of R(X^). The result is: 
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As a technicality, it should be noted that the covariant 
derivatives in Eq. (86) refer to the full 4d space, rather 
than the 3d spatial subspace (as is consistent with the 
fact that R is a 4d scalar density quantity) — despite the 
fact only the spatial components are summed over. 

We now have, in Eq. (86), an expression for the action, 
in terms of R(X^) and its spatial (C 1 ) derivatives, which 
is valid for any set of natural coordinates, X^. Note that 
the time (A) derivatives, d R, do not appear in Eq. (86) — 
as a result of which, the corresponding Euler-Lagrange 
PDE involves spatial derivatives only, and the solution 
R(X fi ) for a given fixed A value is completely independent 
from that for any other A value. The solution R(X^) 
therefore exhibits a constraint in C % but no constraint 
in A — exactly the opposite behavior from that presumed 
in Sec. Ill, and which underlies Postulate 1. Thus — and 
despite being linear (in a covariant sense), and expressible 
solely in terms of the R(X^) field as desired — the new 



PDE is likely meaningless, and in any event, is not at all 
equivalent to Eq. (77). 

A more meaningful result is obtained by introducing 
a slight modification to Eq. (86) — namely, extending the 
summation to include the time component: 
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In natural coordinates, Eq. (87) now gives rise to 
an Euler-Lagrangc PDE with both space and time 
derivatives. The solution R(X 11 ) field is now sensibly 
behaved — at least in most instances. However, even if 
Postulate 1 is now satisfied, Postulate 2 definitely is not; 
the new A derivatives enter into the determination of the 
quantum force components, which thus now depend on 
the future states of the system. Moreover, the PDE is no 
longer second order in A, as there are higher order mixed 
derivative terms that now make an appearance. 

The primary significance of Eq. (87) is that it can be 
interpreted as a completely general covariant expression, 
with X^ an arbitrary set of coordinates. Indeed, com- 
parison with Eq. (85) shows it to be the generally covari- 
ant KG action — apart from an immaterial multiplicative 
constant. This requires that R be identified with <3> (to 
within a multiplicative constant, as discussed), and also 
that the KG constant /z take the value 
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(which is y/2 times larger than the usual assignment). 
The corresponding Euler-Lagrange PDE is the generally 
covariant version of Eq. (84); however, henceforth, we 
will find it convenient to work directly in a X^ = x a 
incrtial coordinate frame. 

The major lesson learned from this exercise is simply 
that the KG PDE fails to satisfy Postulate 2. In accord 
with the discussion in Sec;. IV A, it would appear that we 
can therefore dismiss the KG theory, as being unable to 
provide a viable physical description for individual mas- 
sive particles. In retrospect, various KG difficulties, in- 
cluding those associated with causality, 10 ' 34,36 have been 
known in the single-particle context for many years — 
despite which, it is still used as a basis for many-particle 
QFT. In any event, we will discuss some of those draw- 
backs here, particularly as they relate to the trajectory 
formulation. 

With p(x a ) = Q(x a ) 2 as presumed above, &(x a ) must 
be real valued, and so the KG PDE in question refers to 
electrically neutral particles. This case admits travelling 
wave solutions of the form 



<i>(x a ) 

where h u> 2 — h c 2 k • k 



cos [k • x - cot] , (89) 



2 4 

m c . 



22 



For the $(x a ) solutions of Eq. (89), the wave velocity is 
well known to be supcrluminal. If probability is con- 
served as per Postulate 1, and p{x a ) = $(a; Q ) 2 , this 
in turn implies quantum trajectories that must also be 
superluminal — an untenable situation that should not be 
supported in a relativistic theory of massive particles. 

On the other hand, KG theory adopts a different def- 
inition for the probability density. Specifically, the KG 
flux four- vector is defined as follows: 

j a = ^[$d a $*-$*d a $] (90) 

In Eq. (90) above, '*' refers to complex conjugation — e.g., 
<j>* = $ for the present, neutral particle case. In general, 
the KG j a satisfies the continuity relation, Eq. (28). It 
leads to a scalar probability density, p{x a ) 2 = —j a j a /c 2 , 
and also a spatial probability density, j°/c, both of which 
are (in general) different than <&*<!>. Indeed, for neutral 
particles, the KG j a — 0, in keeping with its interpreta- 
tion as an electrical current flux. Thus, for neutral par- 
ticles, the KG flux adds nothing to our understanding of 
quantum trajectories. 

The KG flux continuity condition becomes non-trivial 
in the context of charged particles, for which $(i a ) is 
complex- valued. In this context, one would like to asso- 
ciate the trajectory-based R(x a ) quantity with |$(x a )|, 
as is done in the non-relativistic case. However, this asso- 
ciation does not lead to a KG scalar or spatial probability 
density that is equal to R 2 = $*$, as discussed. On the 
other hand, the KG forms by themselves do lead to their 
own relativistic quantum trajectory ensemble — defined 
via U a oc j a , and propagating under the influence of a 
quantum potential, Q, that turns out to be proportional 
to the 4d Laplacian (d'Alembertian) of R(x a ). 

As sensible as the KG approach to relativistic quantum 
trajectories — as described above — may seem, it is quite 
problematic. 10 This is true even if one restricts considera- 
tion to the positive energy solutions only — thus avoiding 
Problem (1), as discussed in Sec. I. To begin with, the 
d'Alembertian form of the KG Q clearly violates Pos- 
tulate 2, as discussed. Second, the KG flux four-vector 
j a (and hence U a ) may switch from being timelike to 
spacelike — thus implying quantum trajectories that pass 
through a light cone, which is unphysical. Third, and 
even more egregious: the sign of j° may change, during 
the course of the evolution — thus implying a trajectory 
that turns a corner in time, or a particle that dynami- 
cally changes the sign of its charge. Finally, the inter- 
pretation of the spatial probability density in the two 
theories— i.e., the KG j° vs. the TDSE is funda- 

mentally different, with the KG $ alone insufficient to 
specify a relativistic quantum state (the time derivative 
is also needed). 

It is perhaps remarkable that all of these difficulties 
of the KG approach appear to be remedied though what 
amounts to a fairly simple fix: restrict the derivatives 
used to construct the quantum potential so that these act 
only on the 3d spatial simultaneity submanifolds , rather 



than the full J^d spacetime. Because, in the relativistic 
context, the quantum trajectories themselves define the 
simultaneity submanifolds, this restriction also seems to 
imply that a trajectory-based formulation must be used in 
the relativistic context — or at least, is more fundamental 
than a wavefunction-based approach. Though we have 
not yet succeeded in developing a linear relativistic quan- 
tum wave equation, it remains an open question whether 
or not this is in fact possible. In addition to the reason 
stated above, another important reason why it may not 
be possible has to do with the variation of tj- across a 
simultaneity submanifold, and the fact that relativistic 
quantum trajectories evolve at different rates across the 
ensemble. Adopting a KG quantum trajectory approach, 
it is exactly this effect that causes j°/c ^ <&*<!>, as well 
as the other undesirable properties of the KG j a that 
have been discussed. Alternatively, adopting the present 
approach, the resultant relativistic quantum trajectories 
are evidently well-behaved — but likely at the expense of 
the existence of an equivalent linear wave equation. In 
any case, in the non-relativistic limit — where both the 
wave and trajectory approaches are certainly viable — the 
variability across the ensemble of tj-, and the associated 
trajectory evolution rates, disappears. 



VI. EXAMPLES 

In this section, we present explicit analytical forms for 
solutions of Eqs. (65) and (77) — i.e., relativistic quan- 
tum solution trajectory ensembles, x Q (X M ), expressed in 
natural coordinates, X^ = (cT, C) — for several special 
cases for which such solutions have been obtained. 



A. Quantum inertial motion 

The simplest and most obvious example is that of 
quantum inertial motion. In part to justify some of the 
claims made previously within this document (e.g., in 
Sees. II B, IIC, HIE, VC, and VD), we here work out 
this example in some detail. 

We have claimed that X 11 = x a comprise a quan- 
tum inertial solution, where A M is a set of uniformiz- 
ing natural coordinates, and x a is any inertial coordi- 
nate frame, related to x a via a Lorentz transformation 
[Eq. (5)]. Clearly, g^ v = n^, = S* 3 , and 7 = 1. From 
the second and third lines of Eq. (77), respectively, we 
find Q = f a = 0. From Eq. (75), r r = l, so r = T. 

The equation of motion, as expressed by the first line 
of Eq. (77), then becomes Xj- r = Uf = 0. It remains 
only to verify that this equation is indeed correct, for 
our solution ansatz. From Eq. (5), 

4 = U a =cA a >, (91) 

where A" Q / is the inverse of A" a . The components of 
the inverse Lorentz transform tensor, A°v, are constants 
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that do not depend on the spacetime manifold point, p. 
Therefore, Uj- = 0, as required. 



B. Exponentially decaying probability density 

The case where the probability density decays expo- 
nentially is an interesting example — primarily, because 
it leads to Q =/= 0, but fa = 0. The vanishing quan- 
tum force means that the quantum trajectories experi- 
ence unaccelerated straight-line motion — as in the anal- 
ogous, non-relativistic case. However (and quite unlike 
the non-relativistic case), Q ^ implies a time dilation 
or contraction (in this case, contraction) that is also ex- 
perienced by the relativistic quantum trajectories. This 
example thus embodies a clean separation of the two dis- 
tinct dynamical effects that arise from the relativistic Q, 
as discussed in Sec. V C 

Insofar as the quantum trajectory orbits are concerned, 
these should conform to parallel straight lines, as in the 
case of quantum inertial motion. Without loss of gener- 
ality, therefore, we may take those trajectories to be at 
rest, i.e., 



x l (T,C i ) = S l i C i . 



(92) 



Since T must be orthogonal to the C l , Eq. (92) implies 

t(T,C l ) = t(T). (93) 

Also from Eq. (92), we can take r = t. 

With the x a (T, C l ) ansatz above, we find: 



.9oo = —tt 
9m = 



m = 



(94) 



This also leads to 7=1. We adopt, for the exponentially 
decaying density, the explicit form 



/(C) =exp[-2/c-C], 



(95) 



where n is a constant wavevector quantity, with inverse 
length units. From Eq. (60), the relativistic quantum 
potential is then found to be 



Q = 



2m 



K ■ K 



const 



(96) 



Note that Q < everywhere — implying time contraction, 
rather than dilation, as indicated above. This behav- 
ior is reasonable, considering that exponentially decaying 
probability density is associated with quantum tunneling 
and classically forbidden regions. 

The time contraction effect gives rise, via Eq. (75), to a 
constant value of tj- that is larger than one. This implies 
that the proper time r (also equal to t) progresses more 
quickly than the ensemble proper time, T . In any case, 
we now have an explicit form for Eq. (93), i.e., 



t(T) = cxp 



r. 



(97) 



We are now ready to confirm whether Eqs. (92) and 
(97) satisfy Eq. (65). The first part of this equation of 
motion yields 



dx a 



= c exp 



mc 



K ■ K 



u ) 



(98) 



which is correct, given that U a = (c, 0, 0, 0). The second 
part of this equation, i.e. Uj- — 0, is also true, given that 
f a = because Q — const. 

The relativistic quantum trajectory orbits described 
above are exactly the same as for the corresponding non- 
relativistic application. Indeed, these trajectory orbits 
are also identical to those for quantum inertial motion — 
which may now be regarded as the limiting case of ex- 
ponential decay in which the k 1 — > 0. However, the new 
and interesting feature of the relativistic trajectory-based 
theory is the time contraction constant, tj-. In principle, 
the value of this constant enables one to distinguish be- 
tween different exponential decay applications (including 
the inertial motion case, corresponding to tj — > 1) solely 
from the trajectory ensemble itself — i.e., without direct 
recourse to the probability density function, /(C). Such 
a distinction is impossible in the non-relativistic quantum 
trajectory theory. 

As an interesting exercise, we consider a sinusoidal, 
rather than exponential, form for /(C) — even though 
such a form does not yield a legitimate solution, because 
the probability density /(C) becomes negative. This 
form also gives rise to = and Q = const — but now, Q 
is positive, rather than negative. The relativistic quan- 
tum trajectories are still parallel, unaccelerated straight 
lines, but they experience (identical) time dilation, rather 
than contraction. Interestingly, the resultant R(t, x) field 
is very similar to the KG form of Eq. (89), except that 
the wave velocity is less than the speed of light. 



C. Relativistic Gaussian wavepacket 

In the study of the non-relativistic quantum mechanics 
of free particles, the time evolution of the Id Gaussian 
wavepacket is the canonical example. We define this to 
be a TDSE solution for which the probability density has 
Gaussian form at all times: 



p(t, x) oc exp[— a(t)[x — x c (t)Y 



(99) 



Ideally, an analogous relativistic generalization could be 
derived — giving rise to a nice, well-behaved, nontrivial, 
benchmark analytical solution for the new relativistic 
quantum PDEs of Sec. V. As it happens, though, this is 
not possible. 

To begin with, one is forced to make a choice, based 
on the following simple question: what is the most 
singular characteristic of the non-relativistic Gaussian 
wavepacket? Specifically, is this the fact that the ini- 
tial density is of Gaussian form? Or alternatively, is it 
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that the initial functional form — whatever that happens 
to be — is preserved over time? In the relativistic case, 
we can satisfy one or the other of these two properties, 
but not both simultaneously. 

In this subsection, we mostly presume that the 
first consideration is most important — i.e., that we 
are working with a relativistic wavepacket that is ini- 
tially of Gaussian form (specifically, a coherent state or 
minimum- uncertainty Gaussian). As we shall see, such 
a wavepacket disperses over time — as may be expected, 
based on the well-known analogous non-relativistic be- 
havior. However, in no sense is the Gaussian form re- 
tained over time. This would appear to render ana- 
lytical solution of the relativistic Gaussian wavepacket 
intractible — although a curious result discussed at the 
end of this subsection hints that an analytical approach 
may indeed be possible. In any event, we shall, in the 
meanwhile, obtain the solution numerically. Later, in 
Sec. VI D, we will then ask the other relevant question, 
i.e., what functional form is preserved over time (if any), 
in relativistic quantum mechanics? It turns out that such 
solutions do in fact exist, and moreover, that analytical 
expressions for these can be provided. However, these 
solutions are singular. 

Before addressing either of the above questions di- 
rectly, it is worthwhile to consider exactly why it is that 
the Gaussian form — and only that form — is preserved 
over time in the non-relativistic context. Fortunately, the 
trajectory ensemble approach provides a very straight- 
forward answer. Because of Postulate 1, preservation of 
the probability density form requires that initially uni- 
formly spaced trajectories [i.e., C — xo = x(t = 0, C)\ 
remain uniformly spaced over time. From a dynamical 
perspective, this in turn requires that the initial veloc- 
ity and acceleration fields, x (C) and xq(C), be linear 
in C. The io(C) condition is satisfied automatically if 
^o(x) = ^(t = 0,x) is itself a Gaussian [as opposed to 
just po{x)]. In the overwhelming majority of applica- 
tions, the initial ^o(x) is taken to be a coherent state 
Gaussian, for which xq(C) — const is independent of C 
(i.e., all trajectories are instantaneously moving in uni- 
son). The io(C) condition is satisfied if and only if the 
initial force field is linear in C. For a free particle, only 
quantum forces are present, and so a linear force field 
requires that Qo(C) = Q(t = 0, C) exhibit a quadratic 
dependence on C. The only po(x) functional form that 
can give rise to such a Qo(C) is the Gaussian. 

From Eq. (50), an initially Gaussian po(x) implies a 
Gaussian /(C), because 7 = 1 at t = (since C — x is 
presumed, and 7 = 711 because the application is Id). Of 
course, /(C) is preserved over time — but this property is 
not unique to the Gaussian alone; it is always true, due to 
Postulate 1. Rather, the fact that both sides of Eq. (50) 
maintain a Gaussian form for all t implies something spe- 
cial about 7 — namely, that ~f(t, C) = j(t) is independent 
of C. The value of j(t) indicates the extent of Gaussian 
dispersion or broadening present at time t. Thus, if the 
initial Gaussian is taken to be a coherent state, j(t) > 1 



for all t, with the equality holding only at t = 0. 

Let us now apply a similar analysis to the relativistic 
free particle problem. Without substantive loss of gen- 
erality, we presume an initially stationary coherent state 
Gaussian wavepacket of the following form: 

/(C) oc exp[-aC 2 ], (100) 



with t {C) =t(T=0,C) 
x (C) =x(T=0,C) 
dt^ 
W 
dx 

w 



*T,0(C) 



T=0 



Oo(C) 



mc z 



T=0 



and Q (C) = Q(7"=0,C) = 





C 

exp 




-^(a 2 C 2 



a) 



Again, the fact that /(C) per se remains Gaussian over 
time is not special, as this property would be true for 
any relativistic quantum wavepacket, by virtue of Postu- 
late 2. In what sense, then, should the question of the 
preservation of the Gaussian form be addressed? 

To this end, it is natural to consider the usual spatial 
probability density, j°(t,x)/c. However, for fixed t > 0, 
this function cannot be Gaussian. The reason is that 
the linear (with respect to x) velocity field, x(t,x), that 
this would entail would require the exterior trajectories 
in the ensemble to travel supcrluminally. The one excep- 
tion would be the case of parallel trajectories; but in this 
case, the probability density would be uniform or expo- 
nential, rather than Gaussian (see Sees. VIA and VI B). 
Alternatively, it is more natural in the relativistic context 
to adopt the simultaneity submanifolds — i.e., the con- 
tours of constant T rather than t — as the subspaces over 
which the Gaussian form might be preserved. The rel- 
evant probability density quantity is thus p*{X^). The 
a priori reason why this approach might have a chance 
of being successful is because the simultaneity subman- 
ifolds fan outwards away from the x axis, so that the 
interior trajectories are evaluated at earlier t relative to 
the the exterior trajectories (see Fig. 2). The latter are 
thus afforded more time to spread themselves out. 

The condition that the trajectories be distributed uni- 
formly across a given simultaneity submanifold, labelled 
by 7", is that 7(T, C) = 7(7") be independent of C — in 
exact analogy with the non-relativistic case. In order that 
this condition be preserved over time, one must have, in 
addition to the linear initial conditions of Eq. (100), a 
linear acceleration field, where the acceleration is taken 
with respect to T — i.e., 



<9T 2 



= A a + B a C, 



(101) 



where the vectors A a and B a are C-independcnt con- 
stants. In contrast, the Gaussian form of Eq. (100) pro- 
vides a linear initial acceleration field, but with respect to 
the proper time r, rather than the ensemble time T. The 
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two time derivatives are related to each other by a fac- 
tor of tj-, which depends on Q via Eq. (75), and thus on 
C in nonlinear fashion. The initial Gaussian wavepacket 
therefore does not satisfy Eq. (101), and so the initially 
Gaussian form of p*(X fi ) is not preserved over time. 

Once again, we see in this example a manifestation 
of the principal difference between the non-relativistic 
and relativistic quantum theories — i.e., the distinction 
that exists in the latter case between r and T, both of 
which are dynamically relevant. As discussed in Sec. VC 
and VI B, this essentially corresponds to the distinct dy- 
namical roles played by the quantum force, f a , and the 
quantum potential, Q. This distinction can give rise to 
some quite interesting and competing dynamical effects, 
even for the simple case of the wavepacket dynamics cor- 
responding to the initially Gaussian form of Eq. (100). 
For simplicity, we continue to refer to this example as 
the "relativistic Gaussian wavepacket," even though it 
is understood that the Gaussian form is not preserved 
over time. As stated previously, the relativistic Gaussian 
wavepacket propagation can be computed numerically. 
Accordingly, we now describe the procedure that we have 
used to do so. 

Since the problem is "Id" in the sense that there is only 
one spatial coordinate, the 7 spatial metric "tensor" is 
really just a single number — i.e., 7 = 7 = 7 n . Using this 
fact, together with Eqs. (10), (60), (75), (77), and (100), 
one can express the dynamical PDEs explicitly in terms 
of the T and C partial derivatives of the x a — (ct, x), as 
follows: 



_1_ tc dQ 
m 7 dC 
1 x c dQ 
m 7 dC 
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We have applied Eq. (102) to a specific Gaussian 
wavepacket that has already been considered previously, 
in the non-relativistic quantum trajectory context. 22 For 
this example, the relevant parameters are as follows: 
m = h = 1; a = 1/2. We have also chosen coordi- 
nate ranges that are essentially identical to those of the 
previous non-relativistic calculation — i.e., < T < 10, 
and — 5 < C < 5. The C range encompasses almost all 
of the total probability of the ensemble, whereas the T 
range is sufficiently extensive to incorporate a very sub- 
stantial amount of wavepacket broadening. As for the 
speed of light, c, this is here taken to be an adjustable 
parameter, so as to enable consideration of a range of 
behaviors, from the non-relativistic limit (c=100) to the 
ultra-relativistic limit (c= 1). We will mainly focus on 
the C — o CclSC, RS el specific example that manifests very 
substantial relativistic behavior. 

The PDEs of Eq. (102) were solved using Mathemat- 
ica's NDSolve routine, using the method of lines with 25 



spatial grid points, but without any explicit specification 
of boundary conditions. For c = 3, the entire calculation 
took 2.14 s on a 2.7 GHz Intel Core i7 CPU. Of course, 
we have no analytical solution for comparison; however, 
the numerical solution so obtained satisfies the PDEs for 
both ct and x to an absolute accuracy of 1.4 10~ 6 or bet- 
ter, throughout the coordinate ranges indicated above. 
For the c=100 calculation, the computed quantum tra- 
jectories are virtually indistinguishable from those of the 
previous non-relativistic calculation, and the simultane- 
ity submanifolds are likewise nearly perfectly horizontal 
(constant-^) lines. The relativistic solution therefore does 
indeed reduce to the non-relativistic solution in this limit, 
as predicted in Sec. VD. For the ultra-relativistic case, 
c= 1, the relativistic quantum trajectories very quickly 
begin to fan out, after which they also approach their 
terminal velocities quickly. In certain respects, this lim- 
iting behavior is the same as that of an infinitely narrow 
relativistic Gaussian, as discussed in Sec. VD. Numeri- 
cal solutions for various other c values in the 1 < c < 100 
range have also been computed, but these will not be 
reported on here. 

The c = 3 relativistic quantum trajectories, and asso- 
ciated simultaneity submanifolds, are those presented in 
Fig. 2. This example exhibits very marked relativistic 
effects as manifest, e.g., in the curvature of the simul- 
taneity submanifolds, as well as in the extensive distor- 
tion present in the upper corners of the figure, where the 
trajectories are starting to approach the speed of light. 
For this example, the trajectories remain vertical and 
parallel for a little while at the start of the propagation, 
before fanning outward from each other. On the other 
hand, the simultaneity submanifolds begin curving out- 
ward immediately, as soon as T > 0. This is a clear 
manifestation of the two distinct dynamical effects, as 
discussed above, and in Sec. VC. Specifically, the cur- 
vature of the simultaneously submanifolds is due to the 
C variation of tj — which, due to the latter's dependence 
on Q, is in play even at T=0. In contrast, the fanning 
out of the trajectories is due to /q, which influences the 
trajectory dynamics via the second, rather than the first, 
derivative in T ■ [Note that we have returned to using 
a Q subscript when referring to the quantum force, to 
avoid confusion with the spatial scalar probability den- 
sity, /(C).] 

The center-most, C = trajectory in Fig. 2 — which is 
the one with the greatest probability density — also repre- 
sents a striking demonstration of the dynamical distinc- 
tion described above. Specifically, this trajectory expe- 
riences no quantum forces — and therefore no "accelera- 
tion," in the usual sense involving r. As a consequence, 
this trajectory is straight, rather than curved like all of 
the others in the ensemble. On the other hand, this tra- 
jectory does experience pronounced time dilation — more 
so than any other trajectory, in fact, because the value of 
Q is largest there. Moreover, there is a marked variation 
in the magnitude of the time dilation experienced, over 
the course of this trajectory's evolution. This variation is 
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associated with a corresponding variation in Q(T, C=0), 
that comes about due to wavepacket broadening (see the 
discussion below pertaining to Fig. 5). 

Further dynamical insight into the relativistic Gaus- 
sian wavepacket dynamics may be gained from Fig. 4. 
This is a plot of 7 as a function of C, for various fixed 
values of T, for the c = 3 example. The primary feature 
of this figure is that the 7(C) curves increase in mag- 
nitude with increasing T — thus indicating the expected 
wavepacket broadening, with respect to the simultane- 
ity submanifolds. Note, however, that this increase is 
not uniform across C — i.e., the curves are not horizontal 
lines. If they were, this would imply preservation of the 
Gaussian form over time, which as discussed, does not 
occur (except in the non-relativistic limit). Conversely, 
the curvature or bowing of these curves — which for this 
example, is seen to be rather pronounced — is an indica- 
tion of relativistic quantum dynamical effects. 

Regarding the bowing of the 7(C) curves, a curious 
feature may be observed in Fig. 4, which is that the di- 
rection changes over time. The early curves bow upward, 
whereas at later times, the curves bow downward. This 
can be explained, yet again, as a competition between be- 
tween Q and fg dynamical effects. At early times, all tra- 
jectories are moving (really at rest) in unison; their accel- 
eration has not yet had a chance to manifest as a fanning 
out of velocities. However, even at T=0, the local proper 
time evolves at a higher rate towards the exterior fringes 
of the ensemble (i.e., towards larger \C\). The early si- 
multaneity submanifolds thus curve away from the x axis, 
resulting in larger intervals between nearby trajectories 
that that lie further from C= — and thus, in upward- 
bowing 7(C) curves. Over time, acceleration gives rise 
to trajectory fanout, as discussed. However, as the ve- 
locities of the exterior trajectories reach the order of the 
speed of light, further broadening of the wavepacket is 
relativistically hindered, and so the resultant 7(C) curves 
bow downward, rather than upward. 

The dynamical behavior of the quantum potential, 
Q(T,C), is also worthy of discussion. In Fig. 5 are 
presented curves of Q as a function of C, for various 
fixed values of T (the same values as in Fig. 4), for 
the c = 3 relativistic Gaussian wavepacket. At T= 0, 
Q=Qo(C) exhibits the usual negative quadratic pattern 
[Eq. (100)], well known in the context of non-relativistic 
Gaussian wavepackets. The Qo(C) curve is positive in 
the interior of the ensemble and negative in the ex- 
terior, with the turnover, Qo(CVcf) = 0, occurring at 
C rc f = ±y/l/a — ±y/2. These points thus serve as 
the instantaneous demarcation between "classically al- 
lowed" and "classically forbidden" regions. A remarkable 
feature to emerge from Fig. 5 is that these demarcation 
points are the same for all values of T ■ In other words, 
the boundary between allowed and forbidden regions is 
marked by two special trajectories, C re f = ±V2. We 
refer to these as "reference trajectories," because of the 
fact that Q = everywhere along these trajectories, im- 
plying from Eq. (75) that dr = dT ■ The proper time, r, 
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FIG. 4: Spatial metric tensor, 7 = 711 = jn, for the Id c=3 
relativistic Gaussian wavepacket of Sec. VIC, as represented 
in the natural coordinate frame, X M = (cl~,C). Each curve 
represents 7 as a function of C for fixed T = {0, 1, 2 . . . , 10}, 
with increasing T values corresponding to increasing 7. At 
early times, T>0, the curves bow upward, but at later times 
they bow downward. 



as measured along a reference trajectory, may therefore 
be used as the dynamical ensemble proper time, T, for 
the entire trajectory ensemble. 

It should be emphasized that the existence of refer- 
ence trajectories for a given trajectory ensemble is not 
to be generally expected. We have no general existence 
proof — nor, indeed, has their existence been mathemati- 
cally proven even for the relativistic Gaussian wavepacket 
(although for the c = 3 case explicitly considered here, 
Q(T, ±\/2) = has been established numerically across 
all T, i.e. not just for the eleven T values indicated in 
Fig. 5). That reference trajectories do seem to exist in 
the relativistic Gaussian case is therefore quite special — 
perhaps indicating that in some sense, the Gaussian form 
is preserved, after all. This important property may also 
provide insight that ultimately makes it possible to arrive 
at an analytical solution, although to date such efforts 
have not borne fruit. 

Figure 5 also provides — via Eq. (75) — information per- 
taining to the local time dilation and contraction across 
the ensemble. We see that these effects are greatest at the 
initial time, T=0 — perhaps counterintuitively, given that 
the simultaneity submanifold (i.e., the x axis) is least dis- 
torted there. Over the course of time, wavepacket broad- 
enening leads to a reduction in the magnitude of Q(C) 
across all of the trajectories, C, and thereby, to a reduc- 
tion in the time dilation/contraction effect. As discussed, 
the trajectories in the classically allowed region between 
the reference trajectories experience time dilation (with 
the C = trajectory experiencing this to the greatest 
extent), whereas trajectories in the forbidden region ex- 
perience time contraction. 
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FIG. 5: Quantum potential, Q, for the Id c = 3 relativistic 
Gaussian wavepacket of Sec. VI C, as represented in the natu- 
ral coordinate frame, X M = (cT, C). Each curve represents Q 
as a function of C for fixed T = {0,1,2 ... , 10}, with increas- 
ing T values corresponding to flatter curves. Q>0 and Q<0 
represent, respectively, the classically allowed and forbidden 
regions of spacetime. These are demarcated by two special 
"reference" trajectories, C ro f = ±-\/2, along which Q — for 
all time. 



D. Hyperbolic wavepacket solutions 

In Sec. VIC, we have seen that the relativistic "Gaus- 
sian" wavepacket does not preserve its Gaussian form 
over time, in the sense that p* (T, C) is not Gaussian and 
7(7", C) depends on C, for fixed T 7^ 0. This situation 
begs the question: is there some other choice of /(C), 
for which the analogous properties do hold? We show 
in this subsection that such solutions do indeed exist, 
and moreover, that they can be obtained analytically. It 
should be stated at the outset, however, that these solu- 
tions contain singularities. The solution ensembles also 
contain trajectories that move along a light cone. Since 
every simultaneity submanifold includes at least one such 
trajectory moving at the speed of light, any choice for the 
initial 7" value will lead to initial conditions that neces- 
sarily violate the requirements for a viable solution, as 
laid out in Sec. II. This situation by itself provides suffi- 
cient a priori cause to discard these ensembles; hence in 
that sense, it docs not matter that these solutions develop 
singularities at later or earlier times, because they do not 
technically belong to our theory. Nevertheless, these so- 
lutions arc still worth considering because they are sim- 
ple and analytical, and because they serve as asymptotic 
limits for other solution ensembles that are in fact viable. 
For this reason, it will prove instructive to work through 
them in some detail. 

As per the discussion in Sec. VIC, our starting point 
will be the condition that j(T, C) = 7(7") be indepen- 
dent of C. Actually, it will prove beneficial to first con- 



sider the special case where 

7(7") = const = 1 for all 7". 



(103) 



This corresponds to no wavepacket broadening at all, 
in the sense of the simultaneity submanifolds. How- 
ever, such solutions — should they prove to exist — will 
still manifest broadening in (t, x) space, due to the in- 
creasing curvature of the simultaneity submanifolds with 
increasing T (Sec. VIC). In any case, the fact that 7 = 1 
for all T and C, together with Eq. (60), implies that 
the quantum potential exhibits no T dependence — i.e., 
Q(T,C) = Qo{C) for all T. The same must be true for 

.C 1 - 4 Q =i = fa = - d Q/ dc = ~Qc- 

The fact that the force, /q(C), depends only on C im- 
plies that all trajectories in the solution ensemble undergo 
constant acceleration, in the usual proper time sense. In 
standard relativity theory, the general form for a trajec- 
tory undergoing constant acceleration is well known; 44 ' 45 
such trajectories trace out hyperbolic orbits. Adopting 
the initial conditions from the second and third lines of 
Eq. (100), the solution ensemble must thus be of the fol- 
lowing form: 
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Next, we convert proper time to ensemble time, using 
Eq. (75) and the fact that t = r = T = 0. The result is 
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1)+C. 



It remains only to find a Q(C) such that Eqs. (13) and 
(103) are both satisfied. It can be shown that any such 
(5(C) satisfies 

/ Q (C) 2 = -mc 2 .^(C). (106) 

The trivial solution, /q(C) = 0, gives rise to the uniform 
and exponential cases considered in Sees. VIA and VI B, 
respectively. Otherwise, the most general solution is 



f Q (C) = (107) 
apart from an immaterial shift in C. This leads to 

Q(C) = -mc 2 \og[BC] ; t t {C) = BC, (108) 

with B an arbitrary constant. Note that Q(C) is singular 
at the origin. 

Finally, substitution of Eqs. (107) and (108) into 
Eq. (105) leads to remarkably simple explicit forms for 
the resultant 7 = 1 solutions: 
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(109) 
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The simultaneity submanifolds associated with these so- 
lutions are straight lines — all intersecting at the origin, 
which is therefore a singular point. This situation, and 
Eq. (109) itself, are familiar, in the context of construct- 
ing comoving frames for relativistic particles undergoing 
constant acceleration. As discussed in Sec. I, this leads to 
multiple reoccurences of the same spacetime events (the 
origin in this case), and to simultaneous forward (C > 0) 
and backward (C < 0) time propagation. 

Of course, it is exactly this type of situation that we 
sought to avoid with the present formulation. Again, 
though, we emphasize that the Eq. (109) solution is 
strictly speaking not a part of our theory, as it can be 
dismissed based on initial condition considerations alone. 
More specifically, for the T = simultaneity subman- 
ifold, the initial conditions are not well-defined at the 
origin. In effect, there are two trajectories emanating 
from this point, forming a light cone, and thus moving 
at the speed of light. In any case, T= is ruled out as 
a viable initial submanifold. At other T values, the sit- 
uation is even worse, because Eq. (109) is only defined 
outside of the light cone — thus leaving gaps where there 
is no solution at all. As a result, any also fails to 
serve as a legitimate initial submanifold, in the sense of 
Sec. II. 

Having obtained a solution to the 7 = 1 problem — 
albeit a singular one — we next consider the case where 
7 = 7(7"). One obvious solution of this type can be eas- 
ily obtained from Eq. (109), upon realizing that 500 = 
5oo(C) = —B 2 C 2 is independent of T. This is rele- 
vant when one considers that simultaneous exchange of 
(ct,x) and (cT, C) will flip the roles of g o and gu — 
thus leading to a new dynamical solution undergoing 
constant (and in fact, zero) acceleration, but for which 
#n = 7 = 7(T). The analytical form of the new solution 
is: 

t(T,C) = TcoshLAC] (110) 
x(T,C) = cTsinhL4C] 

Now it is the trajectories of Eq. (110) that are the 
straight lines intersecting at the origin, with the simul- 
taneity submanifolds the associated hyperbolae. Also, 
it is now the region inside rather than outside the light 
cone, that is accessible to the system. In any case, this 
solution has the same sorts of difficulties as Eq. (109), 
and can be similarly dismissed. Since the two solu- 
tions, Eqs. (109) and (110) explore complementary re- 
gions of space, one might consider combining them to- 
gether. Such an approach would "almost" work, from 
the perspective of providing reasonable results through- 
out most of spacetime. However, the origin, and the light 
cone itself, would still be problematic — thus again mak- 
ing it impossible to find a suitable set of initial conditions 
for any choice of T. 

As we have seen, the solution ensembles developed in 
this section are not globally viable. Locally, however, 
they may still serve a useful purpose, particularly as 



asymptotically limiting forms. For example, a relativis- 
tic Gaussian wavepacket in the ultrarelativistic limit of 
small c or large a, is characterized by trajectories that 
are observed to originate from a small region in space, 
and to approach their asymptotic straight line forms 
very quickly. The extreme limit of this behavior would 
thus correspond to an initially "Dirac delta function" 
wavepacket, i.e. to Eq. (110). 



VII. SUMMARY AND CONCLUSIONS 

In a series of previous articles, the trajectory-based 
formulation of non-relativistic quantum mechanics was 
developed in some detail. 20-26 What that work offers, 
in effect, is an alternative — to the usual non-relativistic 
theory based on the wavefunction, and the TDSE 
PDE governing its time evolution. Of course, there 
are other alternate formulations — and interpretations — 
of non-relativistic quantum theory, 6-19 some of which 
have been around for a very long time. Almost all of these 
make use of ^, or of quantities derived directly from it. 
The trajectory-based approach stands somewhat apart 
from these others, because it is formulated independently 
of and of the TDSE. Instead, the starting point is the 
path ensemble x(t, C), from which the solution ensem- 
ble is singled out (via the standard Euler-Lagrange pre- 
scription) as that which extremizes the non-relativistic 
quantum action. This leads to the trajectory-based non- 
relativistic quantum PDE, which happens to be formally 
equivalent to the TDSE. Thus, even though the resultant 
quantum trajectories turn out to be identical to those of 
Bohmian mechanics, the theoretical formulation is com- 
pletely different — as Bohm's approach still requires an 
external field, to serve as the "pilot wave" governing 
the quantum trajectory dynamics. 

As intriguing as a wavefunction-free formulation of 
non-relativistic quantum mechanics may be, the fact that 
the trajectory-based PDE is mathematically equivalent 
to the TDSE — and to all other non-relativistic quantum 
formulations, for that matter — suggests that there may 
be no empirical means available that can distinguish one 
such theory from another. They may be all equally "cor- 
rect," insofar as experiment is concerned. This statement 
is not meant to diminish the value of the various compet- 
ing approaches; each has proven to have its own peculiar 
advantages and disadvantages, vis-a-vis the interpreta- 
tion, computation, and formulation of the non-relativistic 
quantum theory. 

On the other hand, the approach presented in this doc- 
ument represents a very significant departure from all of 
these previous efforts. To the author's best knowledge, 
the present work represents the only (evidently) viable 
relativistic quantum formulation for a single spin-zero 
particle with mass. In this context, the Klein-Gordon 
(KG) theory is the first and best known competing for- 
malism, although its problems as a single-particle theory 
are profound (Sec. VE, and citations therein). More to 
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the present point, however, the KG theory and the present 
theory make different empirical predictions; they are not 
mathematically equivalent, and so it should be possible 
to distinguish between them experimentally. Therefore, 
regardless of numerical convenience or interpretational 
elegance, the issue of which theory is "right" is one that 
can in principle be settled in the laboratory. 

No experiment will confirm all of the predictions of KG 
theory — since these include certain unphysical phenom- 
ena such as scattering probabilities that are greater than 
one. 3 ' 10,34 Likewise, the validity of relativity is beyond re- 
proach, and so no experiment performed in the relativis- 
tic limit will confirm the predictions of non-relativistic 
quantum theory — regardless of the particular formula- 
tion that is adopted. Such experiments may, however, be 
able to validate the present relativistic quantum theory, 
or something similar to it. If this were to happen, then 
apart from the obvious direct value, it could also shed 
light on matters of quantum interpretation and ontology 
that might otherwise linger in the realm of metaphysics 
rather than science. 

This prospect bears some further discussion. Work- 
ing from the assumption of experimental validation of 
the sort posited above, one must consider the possibility 
that there might exist other theoretical formulations, not 
yet devised, that would lead to formally identical predic- 
tions as the present theory — as indeed, we have seen to 
be the case in the non-relativistic limit. For example, one 
might well attempt to relativize the other non-relativistic 
formulations, in analogy with what we have done here for 
the trajectory-based approach. Such a strategy is likely 
to fail, however — because the other formulations rely ex- 
plicitly or implicitly on the TDSE for whose relativis- 
tic generalization is presumably the KG equation. Al- 
though the prospect of a viable single-particle relativistic 
linear wave equation is not entirely ruled out, it seems 
quite unlikely (Sec. VE, and discussion below); in any 
case, no such equation has materialized after more than 
eighty years. Without such an equation, however, it is 
not clear how the other, VlZ-based non-relativistic formu- 
lations are to be relativized, as the linear superposition 
postulate is not satisfied. 3,5 

Under such hypothetical circumstances, the 
present formulation — and any associated ontologi- 
cal interpretation — would become the only viable 
single-particle theory able to explain experimental ob- 
servation, as opposed to one amongst several equivalent 
competing theories. As discussed in Sec. II B, the present 
approach indeed suggests its own natural "many worlds" 
interpretation — wherein each trajectory in the ensemble 
constitutes a distinct world, or universe. This is very 
different from "conventional" many worlds theory, 17-19 
however. For one thing, the trajectory worlds do interact 
with each other (which is, indeed, the source of all quan- 
tum behavior); for another, the worlds do not branch 
over time (new trajectories are not "born," e.g., as a 
result of measurement, but exist for all time). A further 
discussion would be out of place in this document, but 



can be found in our previous work, 22 ' 24,25 and in future 
work specifically addressing the quantum measurement 
ramifications of the trajectory-based approach. 

After the above, admittedly speculative discussion, it 
now behooves us to provide a summary of what has actu- 
ally been achieved here thus far. The key ideas, in rough 
order of logical dependence, are listed below: 

• General covariance 

• Trajectory ensembles 

• Simultaneity submanifolds 

• Postulate 1: probability conservation along trajec- 
tories 

• Postulate 2: dynamical forces do not depend on 
future system states 

• Universality of the quantum potential 

• Principle of least (or extremal) action 

These are the most essential elements, leading to a nearly 
unambiguous theory of relativistic quantum mechanics, 
as presented in this work. 

Of course, the Einstein Equivalence Principle, 33,43 or 
more broadly, the principle of general covariance, per- 
vades everything that we do here — restricting possibil- 
ities, and otherwise guiding towards the relativistically 
proper forms for all quantities. The idea of ensembles 
of trajectories that foliate spacetime is borrowed from 
the non-relativistic theory, but is equally natural in the 
relativistic context, where it is arguably even more im- 
portant. The reason is that it leads directly to the con- 
struction of simultaneity submanifolds for accelerating 
particles — a central idea of this work, crucial to all of the 
subsequent development presented here. 

A rigorous, global definition of simultaneity is prob- 
lematic for classical relativistic theory, because a sin- 
gle trajectory can explore only a small portion of space- 
time. It is not clear that any such viable notion existed 
previously — despite certain clever prior developments in 
this direction, involving intersecting light cones emanat- 
ing from the worldline, and the like. In a trajectory en- 
semble approach, however, the construction of simultane- 
ity submanifolds is extremely natural and unambiguous. 
In this respect, the quantum mechanical generalization 
of relativistic theory actually makes things much easier. 

Postulates 1 and 2 provide important insight into the 
kinematics and dynamics associated with the trajectory 
ensemble time evolution. One of the most surprising con- 
clusions (at least for the author) to come out of this work 
is the idea that the KG equation does not satisfy Postu- 
late 2 — as reasonable as both the equation and the postu- 
late might seem. The latter is, after all, satisfied by rela- 
tivistic and non-relativistic classical mechanics, as well as 
by non-relativistic quantum mechanics; it really should 
be expected to hold in the relativistic quantum context as 
well. A quantum application of the Postulate 2 condition 
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logically requires a global notion of simultaneity. Such a 
structure already exists a priori for the non-relativistic 
case, but for relativistic applications, seems to require a 
trajectory ensemble or related approach. In any case, KG 
theory does not provide a simultaneity structure, leading 
to an evident causality paradox. 

In previous work, 24,25 it has been argued that the 
quantum potential should adhere to a universal form — 
which is quite independent of the particular dynamical 
law (e.g., non-relativistic vs. relativistic). This form 
can be derived without reference to the TDSE — although 
the resultant trajectory-based non-relativistic PDE turns 
out to be equivalent to it, as has been mentioned several 
times. In the relativistic context, for which there is no 
previously existing PDE with which to compare, the a 
priori universality of the Q form obviously plays a central 
role in the development of the formulation. Finally, of 
course, the notion of action extremization — so central to 
classical mechanics, relativity theory and the trajectory- 
based non-relativistic quantum theory — is also applied to 
the present relativistic quantum application, placing all 
of these theories within a single, unified framework. 

It is also to be hoped that we may have helped to 
set the stage for a possible integration with GR, and 
with quantum gravity. Certainly, our use of general 
coordinates even in the SR context, lends itself to this 
purpose — as does the inherent use of density fields. We 
have, in any case, already uncovered some quite intrigu- 
ing similarities — and differences — between the gravita- 
tional and quantum potentials, as discussed in detail in 
Sec. VC. Not least of these is the fact that the value of 
Q itself — and not just its spatial gradient — plays a dy- 
namical role in its own right, as seen explicitly in the 
nontrivial example of Sec. VIC. Also quite relevant (in 
the non-relativistic context too) is the fact that very 
narrow wavepackets experience strongly dispersive quan- 
tum forces, causing very rapid broadening via trajectory 
fanout. Thus, on small length scales — i.e., on the or- 
der of the Compton wavelength, but still larger than the 
Planck length — quantum trajectory effects may prevent 
relativistic singularities from occurring, or may otherwise 
play an important dynamical role, e.g. in inflation. 

There are a number of additional areas for future de- 
velopment that might be considered. Some are straight- 
forward, such as the incorporation of external force fields 
(e.g., the electromagnetic field for charged particles), and 
a recasting of the various density and flux quantities into 
a stress-energy tensor form (together with the develop- 
ment of energy and momentum conservation laws). Less 
obvious will be determining how to generalize the present 
theory for multiple particles of fixed number — i.e., the 
relativistic analog of the many-particle trajectory-based 
non-relativistic quantum theory already in place. Here, 
the biggest and most important challenge will be to 
define simultaneity submanifolds for multiple particles. 
This is not so straightforward, and may well prove to 
be the "Achilles' heel" of the present approach — most 
likely related to the oft-quoted claim that relativistic the- 



ories involving interaction energies comparable to parti- 
cle rest masses must encompass the creation/annihilation 
of particles. 3,5 Of course, one might also attempt to apply 
some of the current ideas in a varying-particle-number, 
QFT context. Given the difficulties of the KG equation, 
one wonders whether it is correct to use this as a basis 
for QFT, as is currently accepted practice. How would 
things change if the present approach were adopted in- 
stead? What impact might this have, e.g., on theoretical 
predictions pertaining to the Higgs boson or pi meson? 

Of course, the basic ideas developed here need not be 
limited to massive spin-zero bosons. It is only natural 
that one might consider revisiting the Dirac field for spin 
1/2 particles, as well as the quantized electromagnetic 
field describing photons, armed with the weapons devel- 
oped in the present arsenal. It is not yet clear what such 
investigations might turn up — i.e., whether they would 
lead to new theoretical predictions, as they seem to have 
done in the present case for spin-zero particles, or to 
predictions in substantive agreement with currently ac- 
cepted theories. In the photon context, of course, the 
simultaneity submanifold concept per se would have to 
be discarded — presumably, to be replaced with an anal- 
ogous family of null hypersurfaces. It is not yet clear 
whether this would even be possible. 

Even for the specific relativistic quantum PDEs devel- 
oped here, there is still much work that remains to be 
done, in terms of characterizing these fully. Can an an- 
alytical solution for the relativistic Gaussian (or other 
nontrival) wavepacket be derived? Are the solution en- 
sembles always well-behaved (subject to the conditions of 
Sec. II), or is it possible that singularities can somehow 
develop? Are the PDEs inherently nonlinear, or can they 
be recast into an equivalent linear wave equation? Under 
what conditions is numerical solution of the PDEs unsta- 
ble, and what sorts of algorithms can be used or should 
be developed? 

One very basic property of particular interest has al- 
ready been established. This is simply that the PDE 
order in the spatial coordinates (C) is twice that in the 
time coordinate (T) — even though these equations are 
perfectly invariant with respect to Lorentz transforma- 
tions. This is not surprising in one sense — i.e., in light of 
the theoretical developments presented in this document 
(notably pertaining to Postulate 2). In the context of 
the history of relativistic wave equations, on the other 
hand, it is highly unusual — if not borderline heretical. 
Traditional treatments virtually always presume that the 
space and time orders must be identical, in order to treat 
these coordinates on an equal footing. The present the- 
ory, in contrast, suggests that there is indeed a funda- 
mental difference between space and time — at least when 
these notions are interpreted in an appropriate manner, 
and linked to a single particle/observer. 

In any case, it should be mentioned that the two-to- 
onc ratio of space-to-time orders that characterizes the 
present relativistic quantum PDEs as described above, 
also characterizes the non-relativistic quantum PDEs — 
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both the trajectory-based PDE, and the TDSE itself. 
This is, in large measure, why the present theory reduces 
so seamlessly to these others in the non-relativistic limit. 
One might also do well to recall that it is precisely the 
lack of such an asymmetry in the KG equation — i.e., the 
fact that this equation is second-order in both space and 
time — that is known to lead to the difficulties described 
in Sec. VE, especially that of defining a suitable prob- 
ability density quantity. These problems were well un- 
derstood in the early days of the quantum theory — but 
what is perhaps quite startling is that the same can be 
said of all of the theoretical tools used here to develop the 
present fix (though of course, the trajectory-based formu- 
lation itself is new). In principle, Einstein or Dirac, or 
one of their contemporaries, could have hit upon this very 



same approach. One (at least the author) has to wonder 
how things might have progressed differently, had this 
idea been born in say, 1932, or even 1952, 8,9 rather than 
in 2012. 
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